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Abstract 

Perturbative and nonperturbative QCD corrections to semileptonic B 
decay rates are calculated and given in analytical and numerical forms. 
We present the analytic expression for the double distribution in hadronic 
mass and electron energy in B — > X u ev e decays through one loop order 
of QCD. We then study the leading twist corrections to that quantity 
within the formalism of HQET and employ it in a proposed method of 
measuring V u b with an accuracy of 10%. Our other point of interest is the 
polarisation of the r lepton in B — ► X c ti7 t processes. Analytic formulae 
are given for the one-loop QCD correction to the longitudinal polarisation. 
The resulting moments of distributions are discussed, incorporating the 
HQET corrections found by Falk et al. The polarisation receives only 
a very small perturbative correction. It can be used in extracting the 
quark masses. Experimental considerations led us to the calculation of 
corrections to the polarisation of r with respect to the momentum of the 
virtual W boson. Analytical results for the 0(1/ ml) terms have been 
combined with numerically evaluated perturbative contributions to give 
the corrected polarisation. 
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1 Introduction 



Semileptonic B decays form one of the most vigorously expanding branch of 
the studies of the Standard Model (SM), the current theory of elementary in- 
teractions jy, ^, [| that has stood the test of some 20 years of experimental 
verification. This model is based on the concept of quantum field theory, de- 
velopped at the time of creation of quantum electrodynamics j|, ||, ||. It has 
thus built on and extended the pcrturbativejRJR ||, |To[ |ll[ as well as nonper- 
turbative f§ ([| ^| ^ |l[|§ methods of field theory. For 

all its success, however, SM is far from completely describing the fundamental 
laws of nature. Even if we refrain from listing the nondiscovery of the Higgs 
boson as a liability, we are nonetheless left with a plethora of questions of the- 
oretical nature. One of them is the origin of the values of the masses entering 
the model as well as the no less enigmatic entries of the Cabibbo-Kobayashi- 
Maskawa mixing matrix. The latter are connected with the topical problem of 
CP violation. Whereas there exist theories that attempt to provide a deeper 
insight on these, like supersymmetry and superstrings, the main obstacle often 
does not seem to lie so much with any lack of ideas as with exact measurement 
and precise theoretical prediction. The latter issue is due to the high level of 
complication characterising SM (not to mention any extension thereof) as well 
as the relatively large value of the strong coupling constant. The strong interac- 
tions cause trouble even in the description of the so-called weak processes due to 
the quantum structure of the theory, making their appearance as higher order 
corrections. It is there that the semileptonic B decays belong. The QCD cor- 
rections bring about a change of the decay rate of up to 20%, so it is important 
to take them into account in any description aspiring to realism. 

This paper is concerned with making some progress towards the determina- 
tion of some of the SM parameters by calculating QCD corrections to semilep- 
tonic B decays, which serve as the field of many tests of the model and are 
currently under heavy experimental investigation. Before reviewing the meth- 
ods and previous results, let us delineate the objectives of these calculations. 

One of them is the determination of the V u b matrix element. The entries 
of the CKM matrix are of great import for the understanding of the physics 
behind the weak processes and theories have been proposed to disentagle the 
observed structure |^3|, [24|, |2q| . If the matrix one day proves not to be 
unitary, which assertion would require precision measurements, it will be an 
undeniable sign of the existence of hitherto unknown species of particles. A 
more tangible problem related to these elements is the CP violation. It is due 
to none other than precise measurements |27], |2^] that we now know that CP is 
directly violated. The value of V u b has always been deducted from semileptonic 
B decays and we propose a way to establish it up to an uncertainty of 10%. 

Our effort went also into finding the corrections to the polarisation of the r 
lepton produced in B decays. Polarisation is a quantity blessed, a least poten- 
tially, with a smaller dependence on parameters compared to total rates. This 
is gorgeously shown in the case of the longitudinal polarisation of the charged 
lepton: even though the rate gets corrected by a fifth of its tree level value, 
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the polarisation is never affected by more than one percent. Such a quantity 
is an excellent candidate for determining the quark masses. These masses have 
been notorious for their haziness but the difference between the beautiful and 
charmed quarks is rather well defined. On the other hand, making predictions 
on absolute rates involves the masses themselves and is rather problematic. We 
will have more to say about that in the discussion of nonperturbative correc- 
tions. 

Since we are unable to perform an exact calculation of the semileptonic decay 
rate (or any, for that matter), we must resort to some approximation or other. 
The zeroth order approximation results in the Born term. This is supplemented 
by corrections which can be divided into perturbative and nonperturbative. The 
perturbative ones are valid as long as the decaying particle can be viewed as 
almost free. Such an assumption is viable in B decay processes because of the 
high momenta transferred and the asymptotic freedom of partons in QCD. As 
mentioned above, the ensuing corrections are substantial and their inclusion is 
a must for an accurate prediction. The technique employed in these calculation 
is that of Feynman diagrams. The accuracy achieved in such a way can be 
estimated by looking at the relative size of some higher order quantities or at 
the various scale prescriptions giving what they claim to be the optimal value of 
the gauge coupling. Unfortunately, the evidence is not uplifting for B decays. 
The total rate gets a two loop correction of a size comparable to that of the one 
loop one, and the scale predicted by BLM scheme is rather low, too. 

At the same time, use has been made of the natural small parameter avail- 
able in the description of the B decays, which is the ratio of the QCD scale to 
the released energy, the latter of the order of the b quark mass. The methods 
exploiting this fact are known as nonperturbative. The early phcnomenological 
approach to explicitly deal with B decays was proposed by Altarelli et al. Q , 
who considered the probability of the b quark inside the meson to have a given 
kinetic energy as a function of this energy. The decay was to be calculated inco- 
herently, which is plausible due to the large energy released during decay. This 
attempt, referred to as ACCMM, was later criticised by the proponents of the 
heavy quark effective theory, HQET, as failing to capture QCD-specific features 
of the decay process. HQET was developped Q as a systematic expansion in 
terms of phenomenological parameters of the order of inverse quark mass. Later 
on, another type of expansion appeared [2^| , with the first nontrivial term corre- 
sponding to the leading twist approximation as used in deep inelastic scattering. 
The latter method emerged after the realisation of the difficulties that the mass 
expansion had had. In the end, the discrepancy between the phenomenological 
method and HQET has been seen diminishing. In our calculations we employ 
the HQET methods, but we take an interesting look at the kinematics of the 
process, using a generalised ACCMM model. Neither method is too illuminating 
about the mass of the quark in question. In HQET, the mass expansion can be 
challenged by questioning the notion of quark mass as a physical quantity. In 
ACCMM, this problem is even less well defined because of the arbitrary charac- 
ter of the incoherent integration. Of course we have to make specific predictions 
so we will state our assumptions clearly in Sec.H. 
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Our work is a continuation of the research done in this field. The semilep- 
tonic decays of B mesons have been investigated for some time using both per- 
turbative and nonperturbative methods. As regards the total decay rate, the 
analogy to muon decay made it possible to use the old results pi, which were 
based on the calculation of radiative corrections performed in H. The 0(a s ) 
corrections to the lepton spectra were first given in a correct, analytical form in 
[ p2| (while numerical estimates were given in p3| ) and then distributions taking 
into account the b quark polarisation were found |34j| . The lepton mass was 
taken into account in |j5| Em| |37[| , allowing to consider semitauonic decays. In 
[ j35l , the one-loop order correction was also found to the longitudinal polarisa- 
tion of the t lepton in the rest frame of the virtual W boson. The polarisation 
turned out to be only slightly affected by the radiative corrections. We have 
found the analytic formulae for the longitudinal polarisation in the rest frame 
of the decaying quark in an analytic form, which is presented in this paper. 

At the same time, the non-perturbative corrections were calculated. As the 
heavy quark effective theory arose together with the application of the operator 
product expansion, inclusive lepton spectra were considered, first in The 
calculation of the average hadronic mass squared discussed in that paper already 
suggested that the expansion in inverse powers of the decaying quark mass may 
start with the quadratic, rather than linear, terms. This statement, which in 
fact needed to be qualified and is now known as Luke's theorem is true for 
the inclusive leptonic spectra, as is seen in particular from the results presented 
here. The actual non-trivial HQET corrections to the leptonic spectra were first 
found in |39| and the work was then extended to double differential spectra in 
electron energy and leptonic invariant mass squared [^C], 41 . All these results 



neglected the mass of the final charged lepton. The latter was considered in 
[ [f2| B B , where the inclusive lepton spectra were calculated. In ]4fj], the 
longitudinal polarisation of the r lepton was found through order 1/mf. We 
apply the inverse mass expansion to evaluating the nonperturbative corrections 
to the polarisation of the tau lepton with respect to the momentum of the 
virtual W boson. The latter is advantageous from the experimental point of 
view [47| . This HQET calculation shows some apparent divergencies, which 
however cancel. A similar, although different, problem was addressed in |48[ . 

The determination of the V u b matrix element is also about finding the per- 
turbative and nonperturbative corrections. We have already mentioned the 
basic papers ||, 0, [52| concerning the perturbative part. On the other hand, 
the HQET provided the inverse mass expansions J3l], 39, [yj. However, it 



was clear from the beginning [ p9| that this kind of expansion suffers from singu- 
larities close to the high endpoint of the lepton energy spectrum. Usually this 
problem may be partially solved by resorting to moments. However, in the case 
of the charmless decays, it is exactly this part of the spectrum that is cleanly 
distinguishable from the much more copious charmed decays and therefore we 
would rather it were treatable on a more detailed basis. One method designed 
to sidestep these difficulties was the so-called leading twist expansion. We take 
this approach and propose a quantity that will enable measuring V u b with an 
accuracy of 10%. 
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This dissertation comprises several results regarding the polarisation in semilcp- 
tonic B decays and the determination of the V u b matrix element. Apart from 
the leading twist convolution of the b — > uei7 e rate, these results have been 
published before. The one-loop QCD correction to the longitudinal r polarisa- 
tion has been published in [[l9|. The paper p3] contains the Born term of the 
polarisation of the r lepton with respect to the momentum of the virtual W 
boson, while the 0(m%) corrections to it have been calculated in plfl . Lastly, 
the double differential distribution in terms of the hadronic mass and electron 
energy has been found through one loop order in f5^| . 

The paper is organised as following. We start with the kinematics of the con- 
sidered processes, in SecJ^. Then we discuss the types of corrections involved in 
the calculation, Sec. ||. At the end of that section, we devote some space to the 
discussion of the Fermi motion inside a B meson. In Sec. |] we present the one 
loop correction to the double distribution in hadronic mass and electron energy 
Sec. H is devoted to the discussion of this quantity in the context of measur- 
ing V U b- It presents results taking into account the leading twist corrections. 
The polarisation of the t lepton in B — > X c ti7 t decays is considered in Secs.|| 
and ^ First, we give the one loop corrections to the longitudinal polarisation 
and then, in Sec. ^, the nonperturbative corrections to the polarisation with 
respect to the momentum of the virtual W boson are given in an analytic form, 
while the inclusion of the corresponding perturbative contribution is performed 
numerically 
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2 Kinematics 



In this section, we wish to exhaust all the kinematic information about the 
processes that are going to be discussed further on. Except for the calculation 
of the 0(l/m%) corrections to the polarisation of the r (see SecJ^), the results 
presented below are found in the rest frame of the decaying b quark. While this 
is rather obvious for the perturbative corrections, the leading twist corrections 
which are used also assume a static quark. The Fermi motion is reflected solely 
in the variable effective mass of the quark. One notes, moreover, that even the 
inverse mass expansion leads ultimately to kinematical ranges identical to the 
partonic ones. 



2.1 Variables 

The four-momenta of the particles are denoted as following: 

Q for the decaying quark, 

q for the final quark (charm or up), 

I for the charged lepton, or 

r specifically for the r lepton, 

v for the antineutrino, and 

G for the final real gluon. 
The system of the final quark and gluon has its momentum denoted as 
P = q + G, while the leptonic momentum is written as W = I + v. Then the 
scaled masses of the final particles are written in terms of the parameters defined 
as 

t 2 q 2 
TOj m£ 

The particles are assumed to be on shell. The following scaled variables have 
been employed in the partonic part of the calculations: 

2^ {l + v? {g + Gf _2v 1 P_ 

m b mf m% m b 

where v stands for the B meson four-velocity in its rest frame. For the pertur- 
bative calculation, the rest frame of the meson and the b quark coincide and 
then one can write, 

_ 2(Q-Z) {l + v? _ (q+Gf 2Q-P 

X 2' 2' Z 2' X P — 2 I J 

ra£ m£ 

where x stands for the charged lepton energy, t for the invariant leptonic mass, 
z for the hadronic mass, and x p for the hadronic energy, all scaled to the units 
of b quark mass. We also define the scaled neutrino energy, x u , but we reserve 
this symbol for the energy in the B meson rest frame. Denoting the meson 
four- velocity as v, we have 

x v = J - . (4 



7 



The charged lepton is described by the light-cone variables: 

r± = ^(x±^-4r,). (5) 

Thus, the system of the c quark and the real gluon is described by the following 
quantities: 



Po = 


\{l-t + z) 


(6) 


p 3 = 


yjp 2 -z = ^[l + t 2 +z 2 -2(t + z + tz)]i, 


(7) 


P±(z) = 


P (z)±P 3 (z), 


(8) 


y P = 


2 P-(z) ^ 


(9) 



where Po(z) and P 3 (z) are the energy and the length of the momentum vector 
of the system in the b quark rest frame, y p (z) is the corresponding rapidity. 
Similarly for the virtual boson W: 

W Q {z) = \{\ + t-z), (10) 

W 3 (z) = yjwg-t = l -[l + t 2 + z 2 -2{t + z + tz)] 1 ' 2 , (11) 

W±(z) = W (z)±W 3 (z), (12) 
. . 1 W+(z) , W+(z) 

^ w -2 ln w^y- ln ^r- (13) 

Kincmatically, the three body decay is a special case of the four body one, with 
the four-momentum of the gluon set to zero, thus resulting in simply replacing 
z = p. The following variables are then useful: 

po = P Q (p) = -{1-t + p) , P3 = P 3 {p) = 'p 2 - p, (14) 
P± = P± (p) = Po ± P3 , w± = W±(p) = 1 -Pzp, (15) 
Y p = y p (p)= 1 -\n P -± , Y w =y w {p) = \\n^±. (16) 

We also express the scalar products in terms of the variables used above, so in 
the units of the b quark mass one gets: 

Q-P=\{i + z-t) , T . v = l{f-n\ 

Q-v= ^(1-z-x + t) , T-P=^(x-t-T)), 

Q- T =^x , v-q=^(l-x-z + ri). (17) 
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2.2 Kinematic boundaries 



The boundaries of the kinematic variables depend on the process which we are 
dealing with. For the partonic part, one can distinguish between three and four 
body final states, the latter including a real gluon. A still different set of limits 
is obtained if one allows for Fermi motion, so that the hadronic entities enter 
the game. This motion is taken into account in the computation of the 0(l/m%) 
and leading twist corrections within the heavy quark effective theory. It is worth 
noting that these expansions have the net effect of formally leaving the b quark 
static. 



2.2.1 Partonic limits 

The scaled electron energy x can vary within the limits 

2\A?<x<l. (18) 
With a fixed value of x, the limits of t are 

V < t < r+(l - — P —) = t 2 . (19) 

However, for the case of a final charm quark this range is conveniently split into 
two regions. In one of them, which we call A, both three-body and four-body 
processes are allowed, while in the other, B, only four-body final states are 
possible. The region A is delimited by 

h=T-(l-— P — )<t<t 2 , (20) 

1 — T_ 

so that the remaining room of the phase space, B, corresponds to 

V<t<h. (21) 
Conversely, if x should vary at a fixed value of t, the boundaries read: 

V < t < (1 - Vp) 2 , w- + — < x < w+ + (22) 

W- w+ 

for region A, and 

r]<t< y/rj{\- 1 P _ ), 207 < x <w_ + ^-. (23) 

for region B. The upper limit of the mass squared of the final hadronic system 
is in both regions given by 

z max = (l-T+)(l-t/T+), (24) 

whereas the lower limit depends on the region: 

p in Region A . . 

(1 - t_)(1 - t/r-) in Region B ^ ' 
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2.2.2 Partonic limits - massless case 

For the b — > uei7 e decay, the masses of the final particles may be treated as 
vanishing. Then the above kinematic limits simplify considerably The electron 
energy x varies within the range of 

0<x<l, (26) 

and for a fixed x, the invariant mass of hadrons can be 

Q<z<l~x. (27) 

If we fix z first, it can take on values from 

< z < 1, (28) 

and then the limits on x are 

0<x<l-z. (29) 

We will also use a triple distribution in terms of electron energy x, hadronic 
energy x p and hadronic invariant mass z. With a fixed value of x, 

l~x<x p <2 — x, max{0 7 x p — 1} < z < (1 — x)(x. p + x — 1). (30) 

As mentioned above, for massless final particles, the distinction between the 
region where only four body final state is allowed and the rest disappears when 
speaking in terms of variables x,t,z fixed in this order. However, with the 
hadronic energy x p used as a variable, an analogous division holds. It is now 
reflected in the fact that the values of x p > 1 correspond to four-body final 
state. 

2.2.3 Meson kinematic boundaries 

In the evaluation of the 0(1/ ml) corrections, it is necessary to work out the 
kinematics in the rest frame of the decaying meson. Let us denote the four- 
velocity of the decaying quark by v. The neutrino energy is no longer fixed 
by the other variables, as it is in the partonic calculation. Rather, it becomes 
another kinematic variable, see ||2j for the boundaries. The partonic value of 
the neutrino energy is given by the relation 

x vQ = l+ t — x — p. (31) 
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3 Evaluation of corrections 



3.1 Perturbative corrections 

Let us denote the final c or u quark generically as q and the charged lepton as 
/. The QCD-correctcd differential rate for the b — > q + /_ + vi reads: 

= dT^ +dTf 3 + dT^ A , (32) 

where 

dr± = G 2 F ml\VcKM\ 2 M± 3 dn 3 (Q;q,T,v)/n 5 (33) 
is the Born approximation, while 

dr± 3 - la s G 2 F ml\V C KM\ 2 M± 3 dll 3 (Q;q,T,is)/ir 6 (34) 

comes from the interference between the virtual gluon and Born amplitudes. 
Then, 

dr f A = \a s G 2 F ml\VcKM\ 2 Mt A dni{Q-,q,G,T,v)/TT 7 (35) 

is due to the real gluon emission, G denoting the gluon four-momentum. Vckm 
is the Cabibbo-Kobayashi-Maskawa matrix element corresponding to the b to c 
or u quark weak transition. The Lorentz invariant n-body phase space is defined 
as 

dTZ n (P;p 1 ,...,p n ) = S^(P-J2Pi)I[ i ^:- ( 36 ) 

The superscript ± refers to the polarisation of the charged lepton when it is 
taken into account. 

The three-body phase space is parametrized by Dalitz variables: 

dTZ 3 (Q; q, r, v) — —dxdt. (37) 

To be specific, consider the unpolarised rates. The evaluation of the virtual 
gluon exchange matrix clement yields: 

^s(r) = - [H q ■ tQ ■ v + H +P Q ■ vQ ■ r + H_q ■ vq • r 

+ ^p(H + + H_)v ■ t + ±p(H + -H_+H L )[r-(Q-q- v)]{Q ■ v) 



where 



-lH L [r-{Q-q-v)]{q-v) 



H = 4(1 - Y pPa /p 3 ) In Xq + (2p /p 3 )[Li 2 (l - 

p+w+ 

W— 

-Li 2 (l -—)- Yp(Y P + 1) + 2(ln VP + Y P )(Y W + Y p )} 



(38) 
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+[2p 3 Y p + (l-p-2t)ln^p\/t + A, (39) 
H ± = ±[l±(l-p)/t]Y p /p a ±jhiy/p, (40) 

H L = ^(l-ln^p) + ^ln^p+^Y pP3 + ^, (41) 

and the superscript "un" stands for the unpolarised rate. After renormalization, 
the virtual correction is ultraviolet convergent. However, the infrared 

divergences are left. They are regularized by a small mass of gluon denoted 
as Xg- In accordance with the Kinoshita-Lee-Nauenberg theorem |l(], |Tl| ,this 
divergence cancels out when the real emission is taken into account. The rate 
from real gluon emission is evaluated by integrating the expression 

i ji« f-\ g i( r ) g 2(r) B 3 (t) 

Mia() ~Jq~W Q-GP-G + W~Gf 1 { ] 

where 

Bi(t) = q- t[Q ■ v{Q ■ G- 1) + G- v - Q ■ vQ ■ G], (43) 
B 2 (t) = q-T[G-v-q-vQ-G + Q-v{q-G-Q-G-2q-Q)] 

+Q-v{Q-rq-G-G-Tq-Q), (44) 
B 3 (t) = Q-v{G-rq-G- P t-P). (45) 

The polarized case requires an appropriate substition of the lepton four- 
momentum r present in the formulae above by a combination of four-momenta 
specific for a sought polarisation. The expressions defining these relations are 
given in the following sections. It will also be seen that this is true both for the 
Born approximation and the corrections affecting the hadronic tensor. 

The four-body phase space is decomposed as follows: 

dK4,(Q;q,G,T,v) = dzdK 3 (Q; P,T,v)dK 2 {P;q,G). (46) 

After employing the Dalitz parametrization of the three body phase space TZ 3 
and integration we arrive at an infrared-divergent expression. 

The method used in these calculations is the same as the one employed in 
the previous ones [37j [52[ [53). The infrared-divergent part is regularized 
by a small gluon mass Xg which enters into the expressions as ln(Ac). When 
the three- and four-body contributions are added the divergent terms cancel 
out and then the limit Xg — ► is performed. This procedure yields well-defined 
double-differential distributions of lepton spectra as described in the following 
sections. 

3.2 Nonperturbative corrections 

The b quark is confined within a B meson so its decay cannot be completely 
described if one sticks with the perturbative, partonic picture where the quark 
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is treated as an almost free particle. While the large momentum transfer ac- 
companying the decay renders it suitable for such an analysis in a theory that 
is asymptotically free, some effects remain and have to be addressed via non- 
perturbative techniques. The HQET methods have established their status as 
the tool to deal with semileptonic B decays and we are using them here. They 
provide a systematic expansion around the limit of infinite quark masses, where 
new symmetries arise, those of flavour and of spin. The expansion comes in two 
kinds. It can be done in terms of consequent orders of the inverse mass of the 
decaying quark, or according to the twist of the operators included. 

3.2.1 Corrections of order l/m| 

Using the operator expansion technique, one can obtain corrections to the decay 
widths of heavy hadrons which effectively lead to new terms in the hadronic 
tensor appearing in the triple differential decay width, 



dx,,dtdx 2n 3 



■ ( 47 ) 



The hadronic tensor W, related to an inclusive decay of a beautiful hadron H b , 

= (2ir) 3 Y,SHPH b ~q~Px) < H b (v,s)\tf\X X X\J c u \H b (v,s) > (48) 
x 

can be expanded in the form 
W^w = -g^uWi + v^VvWi - ie liUa pv a q p W3 + q^Wi + (q^ + v^Ws ■ (49) 

The form factors W n can be determined by using the relation between the tensor 
W and the matrix element of the transition operator 

T^u = -%J d 4 xe-^T[J c J(x)J^O)] , (50) 

which is 

W^ = --Im<H b \T^\H b > . (51) 

The coefficients W n of (^) have all been found elsewhere, see eg. for 
a complete list. Then the distribution ( p7[ ) can be schematically cast in the 
following form: 

(^T 

- — — — = fx8{x v - x v0 ) + f2S'(xu - x u0 ) + fo8"{x v - x v q) , (52) 
dxydtdx 

where 



x v q = 1 + i — p — x (53) 



is the value of the neutrino energy in the parton model kinematics. The triple 
differential distribution must be integrated over the neutrino energy to give 
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meaningful results. The final lepton energy distribution obtained on two sub- 
sequent integrations may be trusted except for the endpoint region where the 
operator product expansion fails. In the present paper we give the double dif- 
ferential distribution so that the lepton energy distribution has to be obtained 
numerically. The calculation does not show any features unfamiliar from the 
cases of the other known polarisations, although apparent divergences pop up. 

3.2.2 Leading twist approximation 

In certain kinematical regions, the inverse mass expansion does not work due 
to the fact that the formally higer terms do not get smaller. Then it may still 
be possible to perform an expansion in terms of increasing twist of the retained 
operators. This technique is taken over from the deep inelastic scattering meth- 
ods. It is applicable when the assumption is valid that the four-momentum of 
the final state hadrons is almost light-like. This is the large momentum that 
has to occur in a twist expansion. 

Under that assumption, the generic quantity T characterizing the B decay 
can be written as 

r(Si) = J dk+fik+JV^^Sm)]. (54) 

In the above formula, Y partonlc is the sought quantity as evaluated within the 
parton model with the b quark mass set to m& + fc+ = m£. The variables 
characterizing the partonic system are generically denoted as £ while those cor- 
responding to the hadronic system by S. The relation between the hadron and 
parton system variables involves kinematics in which the b quark is viewed as 
static but has a variable effective mass m^. Thus the expression has the form 
of a convolution of the partonic rate with the so-called shape function /(fc + ). 
This function is a nonperturbative entity, defined as 

f(k+) =< B\bS(D+ - k+)b\B >, (55) 

with D+ defining the light-cone component of the covariant derivative. 

In the expression ( |54| ) the partonic rate is evaluated for a quark mass equal 
m£, but this replacement is not made everywhere. The prefactor rajj is left 
constant. We will not discuss this prescription, although it seems admittedly ad 
hoc. Let us just note that in the ratio of decay rates, which is what we study, 
this point may be believed not to cause too much trouble. 

3.2.3 Discussion of Fermi motion 

In this paper, what we use are the perturbative corrections and the HQET 
methods of including Fermi motion effects. The latter are worth discussing for 
a while. When it comes to differential spectra, the leading twist approximation 
is the realistic way of using heavy quark theory. In this approximation, the 
dependence of the kinematics on the Fermi motion is rather formal, with the 
quark effectively staying at rest and having a variable mass. We therefore think 
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it an interesting point to make to consider the kinematics of the Fermi motion 
from a general point of view. Even though we abandon the leading twist ap- 
proximation and the kinematic conditions it depends upon, we are still led to 
non-trivial conclusions. 

The products of the decay, the hadronic system and the two leptons, can 
be characterised by three variables, which one may choose to be the hadronic 
invariant mass squared M\ and the energies of each lepton, Ei for the electron 
and E v for the neutrino. The neutrino energy can be measured indirectly using 
energy-momentum conservation. The task of including the Fermi motion con- 
sists in performing a boost from the rest frame of the b quark to that of the 
decaying meson. Then one can sum incoherently over the various configurations 
of the quark in the meson. This probabilistic approach is based on the fact that 
the momentum transfer in the decay is large enough to make an impulse ap- 
proximation valid. We will thus relate the above introduced variables to the 
ones defined in the b quark rest frame. Quantities in this frame will be denoted 
with an asterisk: P* for the partonic four-momentum, E* and E* for the lepton 
energies. Because of the form of the distribution formulae, it is convenient to 
also employ the scaled variables, 

2ET 2E* P 2 , 

xi = , x v = , z = — o. (56) 

m b m b ml 

In the above definition, the mass mb of the b quark will in general vary as it 
is the effective quark mass which defines the kinematics of the process. How- 
ever, we will not vary the factor of m\ multiplying the decay width since this 
factor, reflecting the phase space suppression, is probably compensated for by 
the Coulomb enhancement. As a matter of fact, this mass dependence will drop 
out of the ratio of the charmless to charmed decays. 

In order to specify the relation between the two frames of interest, one needs 
to parametrise the Fermi motion. We do this in a general way, making no 
assumptions on the specific form of the distribution of the quark or the light 
degrees of freedom. 

Working in the rest frame of the B meson, assume that the heavy quark 
has three-momentum k and energy Eb- Forming a bound state, the quark is 

generally off mass shell and its effective mass is then fj, = \J E~% — \k\ 2 . Neither 
are the light degrees of freedom assumed to have the mass of the light quark. In 
fact, using energy and momentum conservation, we write the four-momentum 
K of the spectator as 

K = {M B - E b , -k) = (E„ -k). (57) 

Then we have set up the framework for performing the boost. Starting from 
the system where the b quark stays at rest, the boost needs to endow it with 
momentum k, which defines its parameters to be 

7=- and 0= A (58 ) 
H E b 
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It is instructive to look closely at the transformation formulae for the different 
quantities. The electron energy in the B meson frame can be written as 

E x = \ Xl {M B -E s + \k\ cosOt) (59) 

where Oi is the angle between the directions of the electron and the momentum 
k. Similarly, the neutrino energy is transformed according to the formula 

E v = ^xi{M B -E s + \k\ cos9 v ) (60) 

with Q v denoting the corresponding angle for the neutrino. Note that in each 
of Eqs.( p9| , |60| ) the reference to the Fermi motion is encoded in one rather than 
two variables. These are, respectively for the electron and neutrino, 

= E s — |fc| cos0;, and k„. + = E s — \k\ cos9,. (61) 

It is easy to verify that this structure is due to the fact that the corresponding 
particles are massless. The other case is relevant for the lepton pair, which usu- 
ally has a non- vanishing invariant mass, yielding the following relation between 
the B meson energy, Ew, and the scaled energy in the quark rest frame, 

xw = \{xi + x u) = — — j (62) 



E w — M B xw — E s xw — \k\ cos6 W \l (1 — x w ) — z. (63) 

Clearly, different values of z make different combinations of E s and |fc|cos#vv 
contribute to this expression. In the case of the massless leptons, therefore, we 
can reduce the Fermi motion dependence to a function of the variable kx,+, X 
denoting a generic particle. This justifies the use of the convolution with the 
light-cone shape function. On the contrary, the invariant mass of the lepton 
system cannot be calculated in this way. It also obvious that the double differ- 
ential spectrum of lepton energies also has a more general dependence since it 
involves two different combinations of |fc| and E s . 

Let us now have a look at the energy spectrum of the hadron system. It 
might at first seem that as long as it can be viewed as massless the light-cone 
dominance is in effect. However, one realizes easily that this distribution is 
directly related to that of the W energy by the relation M B = Ew + E B , E B 
standing for the hadronic energy in B meson rest frame. The arguments applied 
to massless leptons do not work here since the energy of the hadrons is the sum 
of the parton and spectator energy. The latter inevitably is disconnected from 
the boost and contributes pure energy rather than the required combination: 

E H =M B x H +E s -x H k +:H . (64) 

For the hadronic invariant mass, the reduction fails to hold if this mass does 
not vanish: 

M\ =E 2 b z + x H E b E s + \p H \ E b \k\ cos6 H + 0(\k\ 2 ). (65) 
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4 Spectrum of hadronic mass and electron en- 
ergy 

4.1 Introduction 

The spectrum of hadronic mass and electron energy is a quantity that can be 
helpful in extracting the V u b matrix element. This section shows the perturba- 
tive corrections to this distribution which are essential due to their numerical 
value of approximately 20%. We will later on take a more comprehensive ap- 
proach at the spectrum of hadron mass and electron energy, but this partonic 
distribution can be useful to let us see clearly what effects come from non- 
perturbative regime and is also evidently a check on the more sophisticated 
calculations. Moreover, while the more realistic calculations rely on modelling 
in one way or another, moments of distributions are believed to be correctly 
given by partonic results due to parton-hadron duality. 

We postpone the discussion of the current status of the \ V u b\ related calcula- 
tions until the next section, where the leading twist approximation is discussed. 
Here, we present the way we evaluate the perturbative corrections in Subsection 



4.2, then show the analytical results, Subsec. 4.2 and us e th em to evaluate a 
few moments of the hadronic mass distribution, Subsec. 4.4. The contents of 
this section has been published in Q . 

4.2 Evaluation 

The tree level result for this spectrum can be read off from the previously given 
formulae, while the one-loop corrections are evaluated following an analogous 
procedure. This calculation, however, required the evaluation of the integral 
over the three-body phase space of the vertex correction as well as the integration 
of the four-body decay rate. 

As regards the virtual gluon contribution, integration over the invariant 
mass of the intermediate W boson gives the desired contribution to the double 
differential distribution, 

dri ' 3 , = 8(z) { -log X G \10x - 25a; 2 + — x 3 + (10 - 24a: + 18x 2 - Ax 3 ) log x 
1210 dxdz 3 3 

-6a; 2 log e + 4a; 3 log el -I- \{2x 3 - 3x 2 )Li 2 (x) + logx [121 

3 18 

-276a: + 195a; 2 - 40a; 3 - logx(30 - 72a; + 54a; 2 - 12a: 3 )] 
121 443 o 128 o 1, . „ o „ o „ 2l 
18 36 27 6 

-4a; 3 log e)}, (66) 

where we have denoted 

x = l-x. (67) 

While we present the decay rate assuming massless final particles, it is not 
possible to eliminate the final quark as well as gluon mass dependence out of 
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the matrix clement alone. This is removed once the integrated real radiation 
term is added. This rate is conveniently split into terms according as they are 
infrared convergent or divergent. Hence we can write, 

3~ conv ~i~ J~i div- (^8) 



12T dxdtdz 
Integration over t gives 

/ Tdiv dt = T d iv, s 5{z) + T div , c 6(z - A G ), (69) 
Jo 

and we obtain the following expressions for the above integral: 

+ log 5(16 - Axz - 14a; + 2a; 2 + 18z + 2z 2 ) + 12.xz + 16a; - x 2 /z 

-7x 2 + -x 3 /z -7z- — z 2 ] 
3 ' 3 J 

+ \ogx[x- 1 {-lQz - I2z 2 ) + x- 2 (2z + 6z 2 ) - 4x- 3 z 2 /3 

4 10 

+22 - 8x/z - I2xz - 36a; + &x 2 / z + 10a; 2 - -x 3 /z H z" 1 - Az] 

3 3 

+ log 2 x(-16 + 4a;z + 14a; - 2a; 2 - 18z - 2z 2 ) + 5" 1 (-20z - 14z 2 ) 
+x~ 2 (2z + 5z 2 ) - ^-x~ 3 z 2 + ^-xz^ 1 - 2xz + 22x - ^x 2 z~ 1 

-9x 2 + ^a; 3 ;^ 1 + 18z+ yz 2 , (70) 

Fdiv,. = -i{9x 2 e 2 logA G + logA G [;o 5 (l-a;)](-30 + 72a;-54a; 2 + 12a; 3 ) 

+log A G (-30x + 57a; 2 - 22a; 3 ) + log A G log —(18a; 2 - 12a; 3 ) 

e 

+ loge(l + loge)(9a; 2 - 6x 3 ) + (-64 + 156a; - 120a; 2 + 28a; 3 ) logx 
+ (15 - 36a; + 27a; 2 - 6a; 3 ) log 2 x ~ 64a; + 37r 2 a; 2 + 100a; 2 - 27r 2 a; 3 
-100a; 3 /3}. (71) 

The infrared finite part has been integrated with the help of FORM. It gives, 
in the same notation, 

/ J~ conv dt — J- conv,s${z) + J~ conv,c@{ z Aq), (72) 
JO 

with 

T conv ,s = logA G (-3a; 2 /2 + a; 3 ) + -!-(173- 360a; + 240a; 2 -44a; 3 )logx 

48 

+ i(3a; 2 - 2a; 3 ) loge + 7^(346a; - 657a; 2 + 294a; 3 ), (73) 
8 96 
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F conV:C = \ogz(x~ 1 (2z-3z 2 ) + x~ 2 z 2 + 9xz + 16x-9x 2 -2z + 2z 2 ) 

+ logx[x _1 (-4.z + 6z 2 ) - 2z 2 x~ 2 + 57/2 - 5xz - 45a; + 25x 2 /2 

-27 z- iz 2 /2] + log 5 log -(-16 + 2xz + 16a;- 2a; 2 - llz + z 2 ) 

+ l/4[x~ 1 (-58z + 6z 2 ) + x~ 2 (+3z + bz 2 ) - 2z 2 x~ 3 - 72xz 

+ 114x - 3x 2 z~l - 61a; 2 + 2x 3 z _1 + 55z - 9z 2 ] . (74) 

The formulae above suffer from infrared and collinear divergences. Thence the 
gluon and final quark masses subsist as regulators in spite of the limit we have 
taken. Of course, both remnant dependences vanish after the integration over 
the hadronic system mass is performed, which involves summing the virtual and 
real contributions. As a matter of fact, the part we have termed as convergent, 
given by Eq. (f72|), may be integrated on its own and then only the logarithm 
of the final quark mass still remains, reflecting the collinear divergence to be 
cancelled against a similar term in the virtual correction. 



4.3 Analytical results 

The first order QCD corrected double differential decay rate can be written in 
the form, 

12T ^ dxdz = Mx)5{z) + ^ [f?6(z) + /f(x, z)6[z - \ 2 G )} , (75) 

with 

1927T J 

where the first term on the right hand side is the Born approximation, given by 

fo(x) = \x 2 {2, - 2x), (77) 



while 



fi{x) = -4/ (x) log 2 A G + i-(120a;- 291a; 2 + 130x 3 )logA G 

lo 



A 



r (-10 + 24x-18x 2 + 4x 3 )logxlog l 



G 

+-(83 - 196a; + 145a; 2 - 32a; 3 ) logx 
6 

+ -^-(249x - 426a; 2 + 155a; 3 ) - |/ (as) U 2 + 3Li 2 (a;)] , (78) 
lo o 



and 



fi(x,z) = -\ -2f (x)logz+ J- (I20x- 291x 2 + 130x 3 ) + ^(10- 24a; 
z 36 3 



+ 18x 2 - 4a; 3 ) logx] + [(10z + 3z 2 )/x - (z + 2z 2 )/x 2 
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+ ^-z 2 /x 3 ] log 4r + \x- l (-mz - 25z 2 ) + ]x- 2 (llz + 25z 2 ) 
3 x A 2 4 

+x~ 3 (-29z 2 /18) + (101/2 - 17xz - Six + 45x 2 /2 
-31z - 3z 2 /2) logx + (32 - 6xz - 30x + 4x 2 + 29z 

+z 2 ) log x log - + log z(21xz + 32a; - 16x 2 - 9z - 5z 2 /3) 
x 

-20x2 + 101x/2 - 97x 2 /4 + 127z/4 + 283z 2 /36. (79) 

The above formula is easily integrated over either of the variables to give the 
single differential distributions in hadronic system mass or charged lepton en- 
ergy. Then expressions confirming previous calculations j3^, Q are found. 
The evident singularity of this distribution disappears after integration over the 
hadronic system mass. That this indeed is so, can be seen by expressing it in 
terms of the following distributions, 

= lim \ hd{z - A) + log A S(z)j , (80) 
l0g ^ = limfe(.-A) + Ilog 2 A^)). (81) 



The substitution of these to Eqs.( |78]j79| ) results in the formal identification, 

e{z-\%)\ = e(z-\%)(~) -5(z)logX%, (82) 



9(z-)&)-f- = 9(z - ) --log* \ G 5(z). (83) 

Upon application of Eqs. ( |S2| , jS3|) to the correction terms, the latter take on the 
following form: 

f[(x) = i(-10 + 24x- 18x 2 + 4x 3 ) log 2 x 

+t:(83 - 196x + 145x 2 - 32x 3 ) logx 
6 

+ -L(249x - 426x 2 + 155x 3 ) - |/ (x) [n 2 + 3Li 2 (x)] , (84) 
18 o 

A c (x,z) = -2/o(z)(^) + (~) [^(I20x- 291x 2 + 130x 3 ) 

~8x + 6x 2 - |x 3 ) logx] + [(102 + 3z 2 )/x - (z + 2z 2 )/x 2 

+ ^z 2 /x 3 l log 4r + \x- l (-mz - 25z 2 ) + tx- 2 (11z + 25z 2 ) 
6 x z 2 4 

+x~ 3 (-29z 2 /18) + (101/2 - 17xz - 81x + 45x 2 /2 
-312 - 3z 2 /2) logx + (32 - 6xz - 30x + 4x 2 + 29z 
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+z 2 )logxlog- + \ogz(21xz + 32x - 16a; 2 - 9z - 5z 2 /3) 
-20xz + 101x/2 - 97x74 + 127z/4 + 283z 2 /36. (85) 

Clearly, the gluon mass does not enter the integrated distribution, defined 

as 

^* *>- I -a<"> + (86 » 

for which we obtain, 

F 1 (x,z) = (ci + c 2 z + c 3 z 2 + c 4 z 3 )log z + c 5 log 2 z + c 6 z + c 7 z 2 + c 8 z 3 + c 9 . (87) 
The coefficients ci to eg are as follows, 

ci = [(120 - 288a; + 216a; 2 - 48a; 3 ) log x + 120a; - 291a; 2 + 130a; 3 ] , 

c 2 = (-30a; + 4a; 2 + 32) log a: + 32a; - 16a; 2 , 

c 3 = 5^ 1 -is- 2 + i(-6a; + 29)logx+(21a;-9)/2, 

c 4 = x -1 - ^x~ 2 + ^.x~ 3 + ^logx - 5/9, 

c 5 = -x 2 /2 + x 3 /3, 
1 , 

-(-102a; + 37a 
(74a; - 33a; 2 )/4, 

c 7 = j (-28a; - 91) logs + \{^- 29) log 2 a; - 10a: -1 log a- "j^ 1 + 

_ 13, 61 145 
x (l ga;+ — ) - — x + — , 

c 8 = — { - 22 log x - 12 log 2 x - 72x~ 1 logx - 162X" 1 
36 1 

+48x" 2 logx + 83x~ 2 - 16x~ 3 logx - 22x~ 3 + 10l}, 
c 9 = _L{6(83 - 196a; + 145a; 2 -32a; 3 ) logx + 12(-10 + 24x 
-18a; 2 + 4a; 3 ) log 2 x - 12tt 2 x 2 + 8tt 2 x 3 + 498a; 

-852a; 2 + 310a; 3 - 36x 2 Li 2 (a;) + 24x 3 Li 2 (a;) } . (88) 
4.4 Moments 

One way of making comparison between the parton model predictions and the 
resonance ridden experimental data is to consider moments of distribution. We 
define those as 



c 6 = -(-102x + 37x 2 + 37)logx + (30x- 4x 2 - 32)log 2 x + 
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While the zeroth moment corresponds to the electron energy distribution itself, 
the singular part of the double distribution leaves no trace on the higher mo- 
ments. In fact, they are then vanishing in the Born approximation. The first 
five moments are expressed in terms of the following functions, 

2a 

M n {x) = —^m n (x), n > 1, (90) 

which are given by the formulae beneath (see also Fig. [l]): 

nil = (-35/144 + 5a; 2 /8-4x 3 /9 + a; 4 /16) logx 

-35a;/144+ 19a; 2 /72- x 3 /48, (91) 

m 2 = (-449/3600 + 7^/24 - 61x 2 /360 - 2a; 3 /45 + 13x7240 
-13a; 5 / 1800 ) logx - 449x/3600 + 533x 2 /1800 - 9x 3 /40 
+109x 4 /1800 - x 5 /lU, (92) 

m 3 = (-103/1800+ 119x/600- 29a; 2 /120+ 19x 3 /180 + x 4 /120 
-3x 5 /200 + a; 6 /600) logx - 103ir/1800 + 697x 2 /3600 
-653a; 3 /2700 + 49a; 4 /360 - 7x 5 /200 + 47x 6 /10800, (93) 

to 4 = logi(-1313/44100+2x/15- 81x 2 /350 + 59x 3 /315 - 5x 4 /84 

-x 5 /175 + 2x 6 /3l5 - 2x 7 /3675) - 1313x/44100 + 5729x 2 /44100 
-lla; 3 /49 + 857a; 4 /4410 - 781ir 5 /8820 + 311a; 6 / 14 700 
-19a; 7 /7350, (94) 

m 5 = logx(-485/28224+ 23a;/245- 211a; 2 /1008+ 151a; 3 /630 

-95a; 4 /672 + 2.t 5 /63 + 29x 6 /5040 - a; 7 /294 + 31a; 8 /141120) 
-485ir/28224 + 1321a; 2 /14400 - 15879ir 3 /78400 
+33599x 4 /141120 - 4523x 5 /28224 + 979x 6 /15680 
-9809a; 7 /705600 + 163x 8 /100800. (95) 



4.5 Summary 

In this section we have presented the results from [p2[ , where they were found 
for the first time. The analytical formulae have been shown for the double 
differential distribution in hadronic mass and electron energy in semileptonic B 
decays. These corrections can be meaningfully used on their own by calculating 
moments of the hadronic mass distributions, which we have presented, too. We 
will extend this study to a more specific method of determining \V u b\ in the next 
section. 
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Figure 1 : Moments of hadron mass distribution m n as a function of the scaled 
electron energy x. The first five moments are shown as indicated by the legend. 
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5 Extraction of \V u b\ from B decays 
5.1 Introduction 

The elements of the Cabibbo-Kobayashi-Maskawa mixing matrix are central 
to the physics of electroweak interactions. A precise determination of these is 
essential in establishing the unitarity of this matrix and so the completeness 
of the Standard Model, as well as in unveiling the nature of the CP violating 
interactions that remain puzzling. The beauty to up transition element, V u b, is 
one of the smallest and least known of the matrix. The current PDG estimate 
ranges from \V u b\ — 0.002 to \V u b\ — 0.005. It is therefore currently a significant 
problem of particle physics to obtain a better estimate of it. As B factories are 
now under operation, the experimental data can be expected to offer greater 
accuracy so that theoretical prediction should also aim at a higher level of 
sophistication. 

The V u b matrix element has long been extracted from semileptonic B decays 
owing to their relative cleanness. They are basically due to the weak interac- 
tions, but are affected by chromodynamical quantum corrections. The inclusion 
of the latter is vital due to their significant size of up to 20%. The process at 
hand has been studied using perturbative and nonperturbative techniques to 
take QCD effects into account. In our present case, QED calculations are of use 
due to the topology of the diagrams involved. While the total decay rate 
of b to u can be inferred from an analogous analysis of the muon decay, due 
to Berman |8| , numerical calculations of spectra appeared in |53| , and the first 
correct analytical one- loop results were published in p2[ . 

All of the calculations mentioned so far dealt with the perturbative correc- 
tions, which are valid due to the large momentum transfer occuring in the B 
decays. On the other hand, nonperturbative corrections are not negligible, but 
due to their relative smallness can be evaluated using an expansion in the scale 
provided by the ratio of the QCD scale to the transferred energy. A first, phe- 
nomenological approach was proposed in [ [30"| , where the b quark inside a meson 
was assigned a probability distribution describing its Fermi motion. The decay 
was then calculated on an incoherent basis. This method, referred to as AC- 
CMM after the names of the article's authors, has later been dismissed by the 
proponents of a more systematic formulation, the heavy quark effective theory 
(HQET), as falling short of describing some aspects of the decays, allegedly 
accessible with HQET. With a systematic expansion proposed in pi]| , the new 
method was used to calculate the first nontrivial corrections ||| to lepton spec- 
tra. These corrections were expressed as an expansion in inverse powers of the 
heavy quark mass and involved terms viewed as breaking the heavy quark and 
spin symmetry, two symmetries of the effective lagrangian. Double differential 
spectra were considered in jy], while the mass of the final lepton began to 
be included with the papers of Balk et al. and Koyrakh j|3| . 

Aside from these inverse mass power corrections, another technique was 
developped |22[ consisting in a leading twist approximation as familiar from 
deep inelastic scattering. This new approach stemmed from the observation 
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that the operator product expansion (OPE), forming the basis of the power 
corrections, fails at the high endpoint of the electron energy spectrum. This 
regime turned out to require a separate treatment with new nonperturbative 
input which was introduced in the form of the so-called shape function. This 
function describes the light-cone distribution of the b quark within B meson. 
The light cone approximation holds as long as the final hadrons can be treated 
as nearly massless. It is worth noting that at this point the effective theory made 
contact with the phenomenological approach. Indeed, analyses have been done 
| |56| to establish the connexion between the two. At any rate, these essentially 
model-like descriptions appear to be necessary whenever one focusses on spectra 
rather than restricting oneself to moments, which are believed to be reproduced 
correctly from the mass expansion. The shape function does not depend on 
the final states but rather is inherent to the B meson itself. This in particular 
means that information gathered from other processes can be used. Along these 
lines, the case of semileptonic B decays is related |57j] to radiative b — > sj 
decays. Formally, the leading twist expansion corresponds to resumming the 
most singular terms in the inverse mass series so that it is perfectly possible to 
calculate both kinds of corrections and disentangle the different terms in order 
to avoid double-counting any part. One may think of those two methods as 
summing the terms appearing in the full result along different paths, and we 
know which they are. 

In the following section we explain our choice of the method of extracting 
| V u b | . Then Sec. 5.3 is devoted to the articulation of both corrections in the 



desired approximation. Subsequently, we discuss the shape function character- 
ising the nonperturbative input used in our calculations in Sec. 5.4 and then go 



over to the results in Sec. 5.5 



5.2 Method of measuring V u b. 

In our paper, we propose to use the semileptonic B decays in the determination 
of the V u b element. One of the difficulties to be negotiated is the prevalence of 
charmed final states, which exceed the charmless ones by a factor of 100. We 
try to address this problem by imposing an upper cut on the invariant mass of 
the final hadrons. Since the lowest charmed state has a mass squared of ap- 
proximately 3.5 GeV 2 , this value suggests the placement of our cut. Whatever 
the exact value chosen, this kind of cut has consequences on the methods we 
have to apply. As regards perturbative QCD corrections, we use the one-loop 
level spectra, which is the best result available. The question is what nonper- 
turbative effects should be included. As mentioned above, there are at least 
two approaches to the problem. It turns out that the cut on the hadronic mass 
makes the resulting distributions sensitive to the Fermi motion of the b quark 
within the decaying meson. Due to this fact, we include the leading twist cor- 
rected terms. As explained in the section discussing the HQET expansions, this 
effectively consists in performing a convolution with the shape function. Thus 
both the tree level and one-loop terms of the partonic distribution are weighed 
with this function with the mass of the b quark effectively varied. 
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The convolution needed to obtain the leading twist approximation requires 
a choice of the shape function. Since it is of a nonperturbative nature, we need 
to decide on a model with parameters to be fitted from experimental data. As 
regards our present knowledge about the shape function, it is usually expressed 
in terms of its moments and the difference between the meson and quark mass 
A = niB — m b . These are not entirely independent quantities, however. The 
best known moment is the upper limit of the support of the shape function, 
related to A. The first moment vanishes by equation of motion, while the second 
corresponds to the kinetic energy of the b quark in the meson and is proportional 
to the parameter Ai, which is rather poorly known to date. The third moment is 
considered in the literature, but we ignore it altogether as little is known about 
it quantitatively. Summarizing the properties of the shape function f(m b ), we 
can write the following conditions, of which the first is normalisation. 

M B 

f(m* b )dm* b = 1, (96) 

-oo 
M B 

f{m* b ){m* b - m b )dm* b = 0, (97) 

-OO 

M B 2 

f{m* b )(m* b - m b ) 2 dml = --Ai, (98) 



M B 



f(ml)(ml~m b ) 2 dml = --p 1 . (99) 

In these formulae, denotes the effective mass of the b quark. Certainly, as 
one gets to know the moments of the shape function better, it will be of use to 
include them all. This is technically reduced to examining the shape function 
to be used in the convolution for an extended range of parameters. 

For our analysis, we have chosen to use the form of the shape function 
suggested in the BaBar physics book Q, which reads 

/(m*)=/(|) = 7(^^), (100) 

with 

f(x) = N6(l - x)e<*- 1 \\ - x)^. (101) 

We have scanned the results obtained for a couple of parameter sets correspond- 
ing to realistic shape functions. The error estimate of 10% which we have arrived 
at is based upon the range within which the result falls subject to the choice of 
the nonperturbative and perturbative parameters. 

In making our choice of the quantity to be studied, we were also guided 
by consideration of the mass definiton. The quark masses affect the decay 
widths rather substantially, entering with a fifth power in the overall prefactor. 
In the prescription of effective mass convolution we have skipped the point of 
this prefactor. Although it is by no means a proved statement, an analogy to 
muon capture suggests that these factors, determining the total rate, should 
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be left constant and that is what we do. Note that this is not the case in the 
original ACCMM model. Whatever the consequences of the treatment of mass 
be, however, we hope that we can ease this dependence by considering the ratio 
of charmless to charmed events. 

As already mentioned, our analysis involves a few parameters and two cuts. 
The most important cut that is freely adjustable is that on the invariant hadronic 
mass. As will be seen from the results, the greater the cut the weaker the sen- 
sitivity to the shape function parameters gets. On the other hand, the very 
reason for placing this cut, the removal of charmed states, requires it to fall 
in the vicinity of 3.5 GeV. These two arguments are clearly contradictory and 
we would not feel comfortable picking one value for all the subsequent analysis. 
Rather, we have decided to consider this cut as a variable and study the above- 
mentioned ratio for different values of the cut. It will be seen that too low a cut 
would mar the result with wild dependence on fine details of the nonperturba- 
tive shape function, which we do not know that well. Eventually, the range of 
2 GeV 2 through 4 GeV 2 will be studied. 



5.3 Combining perturbative and nonperturbative correc- 
tions 

In order to provide means of extracting the ratio of Vckm elements, one has to 
define a quantity that is to be matched to experimental results. We examine 
the ratio 



r u (M* ut ,E cut ) _ |K 



2 



U ~ T c {E cut ) (W2) 



where the rates arc defined as follows: 



T u {M 2 cut ,E cut )= / AM 2 X / dEi 

JO J E cut 



AM 2 x AEi 



(103) 



AT(B -> X c ei7 e 



Tc(Ecut) = J E m ~ ^ AE^^' ' (1 ° 4) 

In the above formulae, E cut denotes the minimal electron energy, while M 2 ut the 
maximal hadronic invariant mass squared. For the charmed decay, the maximal 
electron energy is given by 

E max = <_Z^ i (105) 

as it is treated in the perturbative approximation to order 0(a s ). Thus the 
expression for the decay to the charm quark is in fact taken as 

{ x ~ p AT 

T c (E cut ) = / — Ax, (106) 
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which involves partonic distributions only. The mass of the b quark which is 
taken for this formula is related to the parameter A describing the shape function 
using the approximation 

m b = M B - A, (107) 

valid up to 0(1/ m 2 ) corrections. Thus the dependence on the shape function 
actually enters both decay rates defining the ratio (102). The rate of the decay to 
X u states is computed taking into account both perturbative corrections to order 
a s and the leading twist non-pertrubative effects due to the Fermi motion, which 
are expressed in terms of the shape function /(m*). The resulting distribution 
is 



dMjm = I dx L dxp L n dz L dm * f{ml ^-a- z x 

8(M 2 X - m* 2 z - x p A*m* - A* 2 )S{E l - ^). (108) 
The partonic rate can be split into the tree level term and the 0(a s ) correction, 

^j^parton ^^parton.O _|_ ^^£_^ppar£on,l (109) 

37T 

and then the convoluted rate can be divided up accordingly. Using the delta 
functions, the integrations involved in these rates can be simplified, yielding, 



dEi / AMI , " 

Ecut Jo dM 2 x m 

Mb/2 px 2 / 2E[ \ (^pP ar * on '^ 



JE cut J2EJMB Jy min X \ X 



dxdy 



(110) 



y mm = WO, l-Mi % - (Af B -2y 

y X {2E l /x)(M B -2Ei/x) S K 1 

for the tree level contribution, while the a s correction contributes 

Mb/2 /*1 p2—x fZmax 9 / 9 rp \ r ]j^parton,l 

dE t / dx dx p dz-f(^) ^ , (112) 

E cut J2E l /M B Ji-x Jz min x \ x J dxdx p dz 



where 



and 



1 ^ x v<\^ (113) 
1 for x p > 1, y J 



^M^-^Ms-mM -jMs^E./x) 2 ^ (i _ x){xp + x _ 



(114) 
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The integral over z in Eq. ( |112| ) has been performed analytically so that we use 
the prime function F(x, x p , z) defined as 



F{x, Xp 7 z) 



dz' 



dr pa 



Tq dxdxpdz' 



F = 



F^z) - Fi(0) + F 2 {z) + F 3 , x p < 1, 

F 1 (z)-F 1 ( ) + F 2 (z) - F 2 (z min ) x p > 1, 



where 

and 

Fi ■■ 



(115) 
(116) 
(117) 



\ ln ^ ~ 4(x ~ 1 + V 2 )' A'] +8(x-l + x p /2) 2 /t 2 

x (36 - 3§xp + 12x 2 p - 3x1/4) + - In iil^ (as - 1 + x p /2)[48 - 78x p 
+45^ - 21x 3 p /2 + 3x 4 p /A + (x - 1 + x p /2)(-42 - 15x p + 27a; 2 , 



-21/8^)] + {Ait + A 2 ) ln ^ + A 3 log(l - t) + A 4 . 

^2 = 12[21n 2 z- (81nx p -7)ln z]v(x — x p ), 
F 3 = (x p ~ x) {2AhiXp(~\ + hv — AvhiXp) 
+v[-16tt 2 - 60 - 48Li 2 (l - x p )]} . 

In the above formulae. 



(118) 
(119) 

(120) 



Ai = 



A 2 



(x p — x)(xpZx/2 + 2\x v xj2 — 6x p v — I2x 2 x + 3x 2 v + 7x 3 x/4) 
+ (2x p zv - 9x p z/2 - 9/10a;pZ 2 - x 2 zv + 7x 2 p z + 12x 2 v + lhx 2 /2 
-21x 3 p z/20 + l7x 3 p v/2 - 105/4^ - 5/4a£u + 53/4^ 
-219/160^), (121) 
(x p - S)(-36a;p(l - x p x) + &x 2 v + 4x 3 x + 122 - 96uln2 + 48w) 
-24a;p + 48a: 2 + 16x 3 p v - 60x 3 p - 2x A p v + 26a;* - 12/5x 5 p , (122) 
2A 2 + 96{x p - x)(l + x-x p ) ln(l + t), (123) 
(yX ip x J y 5 x p z x 24zx + I2zv) + 32x p zv — 81a;pZ — 21/5a;j,z 
-Ax 2 p zv + A0x 2 p z - 81/20a;pz + 36zu + 12z- 6z 2 v + 24z 2 , (124) 



and 



x = 1 — x, 



v = 2 — x — x 



V 



t=Jl- 4z/x*. 



(125) 



5.4 Shape function parameters 



The decay rate distributions depend on the strong coupling constant, a s , as well 
as the parameters characterizing the shape function. The latter are identified 
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with the meson-quark mass difference A and the quark kinetic energy inside the 
meson, — 3A1. These physicai quantities make their appearance in the shape 
function according to the BaBar ansatz p8| ], 

/(m*) = ) ee (126) 

with 

f(x) = N6(l - x)e<*-V(\ - x)^-^. (127) 
The normalizing factor N is given by 

N= fm' (128) 

and the parameter c is related to the physical ones by 

-3A 2 



A, 



(129) 



We have considered a set of parameters A, Ai. They are believed to fall in the 
ranges 

0.3 GeV < A < 0.5 GeV, -0.6 GeV 2 < Ai < -0.1 GeV 2 . (130) 

However, we have also imposed the condition that the mean kinetic energy of 
the quark should not exceed considerably the meson-quark mass difference. For 
a s we have used 

a s e {0.2,0.3}, (131) 

For each value of a 8 , the following sets of non-perturbative parameters have 
been examined: 

(132) 
(133) 
(134) 
(135) 
(136) 
(137) 

The parameters of the shape function as well as the value of the coupling con- 
stant have a twofold effect on the examined ratio r since both the numerator 
and the denominator of the expression depend on those numbers. The decay 
rate to charmless states has the dependence encoded in the convolution with 
the shape function and the perturbative corrections. The charmed rate depends 
on the A parameter through the ratio of b to c quark masses. The Born rate 
normalised to To of b — > c decay is 

T = 0.44064 for A = 0.3 GeV, , (138) 
T = 0.49938 for A = 0.5 GeV, (139) 



A 


= 0.3 GeV, 


Ai = 


-0.1 GeV 2 , 
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= 0.3 GeV, 


Ai = 


-0.2 GeV 2 , 
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= 0.3 GeV, 
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-0.25 GeV 2 
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where the electron energy cut was set to E cu t = 0.5 GeV. 

The shape function corresponding to the 6 sets of parameters is plotted as a 
function of the effective b quark mass in Fig.||. For clearer distinction see Figs. 
| and |. 

5.5 Results 

5.5.1 General result on \V u b\ accuracy. 

If we choose M\ = 4 GeV 2 as the hadronic mass squared cut, then the un- 
certainty of the extracted value of \V u b\ will be about 10%. The source of this 
uncertainty lies with the parameters that go into it. The quantity r which we 
have chosen for use in the determination of the ratio of IKjb/V'chI depends on 
the strong coupling constant, the quark masses and the parameters A and Ai 
entering the shape function. This makes for quite an involved structure of the 
various contributions. Before we specify and discuss these, let us briefly recall 
the range of the survey that has been performed. In order to make a prediction, 



we have examined the various sets of parameters as defined by Eq. ( 132 - 137 ). 
The six sets of nonperturbative parameters span a conservative range of the 
estimates for the meson-quark mass difference (A) and the kinetic energy of the 
b quark Fermi motion (Xi). In addition, the strong coupling constant was taken 
to be 0.2 or 0.3 for each set, while the electron energy cut was kept constant 
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Shape function 




effective mass of b quark 



Figure 3: Shape function for A = 0.3 GeV. 



at all times and set to 0.5 GeV. This particular cut is arguably of minor im- 
portance to the analysis and has been made to provide for the identification of 
the semileptonic mode. More vital is the remaining upper cut on the hadronic 
invariant mass. We look closer at the range of 2 GeV 2 to 4 GeV 2 of this cut. 

By examining the variation of the observable r due to the changes of the 
parameters we can obtain an estimate of the accuracy that can be achieved. 
The form of the function r(M^) obtained for all the studied sets has been 
presented in Fig. |ll[ The value of r(M 2 ut = 4 GeV) can be approximately 
written as 

P = 1.74 + 0.3^ + 0.2^ + O.OS^i. (140) 
a s A Ai 

These numbers give a flavour of the dependence on the parameters. Let us now 
turn to the identification of the effects brought about by the nonperturbative 
and perturbative corrections. Since it is the hadronic invariant mass cut that 
is variable in our analysis, we will look at the value of the examined ratio r as 
a function of this cut, r(M] c ). We will look at the partonic result at the Born 
level and at the corrections it receives from the perturbation theory as well as 
HQET leading twist approximation. 
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Figure 4: Shape function for A = 0.5 GeV. 



5.5.2 Partonic results. 

Although it might appear that partonic results, understood to treat the decay 
as that of a free quark into a free quark and a pair of free leptons, have little 
in common with nonperturbative data, the latter intrude via the quark mass 
ratio. This is a simple consequence of the difficulty that arises as soon as one 
tries to treat quarks as free, ignoring confinement. Thus, while the 5->h decay 
rate is truly devoid of any reference to nonperturbative quantities in this zeroth 
approximation, the decay to charmed states will bear a mark of the A parameter, 
which sets the b quark mass. The values of the scaled decay rates have been 



quoted in Eqs.(138,139). With the mass of the b quark fixed, the ratio can also 
be thought of as fixed due to the fairly well known mass difference between 
b and c quarks of 3.4 GeV. Curiously enough, we must conclude that if we 
only considered the partonic result, the fact that we are calculating ratio rather 
that the rate itself would be counterproductive to precision. But of course such 
reasoning would be unrealistic as the corrections are essential to accuracy and 
must be included. We do not present any plots for the partonic quantity of 
r(Mj c ). Since the invariant mass of the hadrons is always zero in this approach, 
such plots would have to be series of horizontal lines. However, beneath we give 
the parton model values of r obtained assuming A = 0.3 GeV and A = 0.5 GeV. 
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Figure 5: B —* X„ decay rate at tree level. 
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Figure 6: B — > 1 M decay rate at tree level in the whole range of hadronic mass. 
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Partonic ratio r as a function of A 
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Figure 7: The partonic value of the ratio r of charmless and charmed decays as 
dependent on the A parameter. 



They are clearly the reciprocals of the scaled rates stated in Eqs.(13£. 139| ). 



r = 2.2694 for A = 0.3 GeV, 
r = 2.0025 for A = 0.5 GeV. 



(141) 
(142) 



The dependence of the partonic value of r on the parameter A is nearly linear 
in any reasonable interval, as seen in Fig. 0. 



5.5.3 Convoluted partonic results. 

The rather dull hadronic mass spectrum delivered by the partonic model gains 
colour when the Fermi motion is accounted for. The distribution gets indented 
at the lower end of the spectrum because now the hadronic mass is taken to 
incorporate also the 'brown muck' or the light degrees of freedom, or, couched 
in the language of phenomenological models, the spectator. This kinematical 
argument is relevant even though we are dealing with a formally developped 
approximation to QCD, which renounces model dependence. This is because 
the kinematics of this approximation is clear. The b quark is static with its 
mass, as it were, bulging or shrinking. The conservation of four-momentum 
then leads to a unique expression for the hadronic invariant mass, 

M\ = m* 2 z-x p A*m* b - A* 2 . (143) 
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Ratios at tree eve 
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Figure 8: Ratios charmless/charmed at tree level. 



This formula implies that the typical scale on which the effects of including the 
Fermi motion are expected to be noticeable is set by Am;,, which is around 2 GeV 
and whose variations are discernible. The extent to which the light constituents 
of the B meson exert influence over the hadronic mass distribution is gauged by 
the parameters of the shape function. In particular, the departure of the b quark 
mass from the mass of the meson, as indicated by A, shows positive correlation 
with the shift of the distribution, see Fig. || showing the rates and Fig. || for the 
ratio r(M^). Note that the probabilistic nature of the convolution implies that 
barring the effects of the imposed cuts one should recover the partonic results 
upon phase space integration. That this does happen is visible from FigJ| where 
the different sets tend to a common value as the hadronic mass is cut sufficiently 
high. 

5.5.4 Convoluted one-loop corrections. 

As explained in Sec. [5.3] , the perturbative corrections are combined with the 
nonperturbative ones in essentially the same way as the tree level partonic term. 
It is well known that one-loop corrections themselves yield a singular hadronic 
mass distribution, characterised by a Dirac delta peak at zero mass and a con- 
tinuous but divergent part. The Kinoshita-Lee-Nauenberg theorem [|l0| 
ensures that the divergences cancel once they are properly regularised. In a 
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strict perturbation theory, final results are only obtained upon integration over 
some final range of the invariant mass, enclosing an interval around zero. The 
convolution prescribed in the leading twist approximation performs this inte- 
gration so that no infinity enters the one-loop corrected and convoluted terms, 
as seen in Fig. ^. However, the large negative correction remains in the convo- 
luted spectra as a dip close to the lower end. As with the tree level result, the 
corrections are affected by the Fermi motion in a degree that depends on the 
size of the applied hadronic mass cut. Clearly, the low end of the spectrum is 
very sensitive to the details of the shape function, while this influence is washed 
out above the rough scale of Amj,. Again, we recover the partonic result after 
inclusion of all phase space, see Fig. 

5.5.5 Perturbative and nonperturbative corrections. 

Summarizing, it is seen that although low cuts on the hadronic mass are favoured 
by the requirement of cutting off the charmed admixture, they leave the resulting 
distribution uncomfortably dependent on the Fermi motion. This motion is 
known via the moments of the shape function and since this knowledge is not 
very accurate we cannot rely on the details of it. Therefore we suggest that 
to be conclusive one must take the cut at least at the level of M\ = 2 GeV. 
On the other hand, to prevent the charmed states from excessive proliferation 
and swamping the charmless events, we do not consider values higher than 
M\ = 4 GeV. 

Deciding upon these ranges, we can come to conclusions about the attainable 
accuracy of \V U b\ as measured in the way we propose. To this end, we have 
collected the calculated values of the ratio r(M^) in Fig. |ll|. The spread of the 
ratio for a given cut reflects the uncertainty ascribed to the measurement of the 
square of the matrix element. Thus, if one sets the cut to 4 GeV, one can hope 
of extracting \V u b\ with an uncertainty of around 10%. 

5.5.6 Discussion of uncertainty. 

The above conclusion regarding the precision we assign to our method aims to 
be conservative and it is based on a scan of a wide range of parameters entering 
the distributions. However, it is possible to improve on the accuracy of these 
parameters. It might be useful to exploit the analogy to the b — > s"f decays and 
extract the nonperturbative information from them, as suggested in J57j. On 
the other hand, the question remains whether one can use other observables of 
the semileptonic decays to help establish parameters. This latter method has 
the advantage of eliminating possible systematic errors due to different kinds of 
considered processes. 

We have taken a look at the differential distribution of the hadronic mass. 
This is to be distinguished from the integrated distribution discussed so far. As 
already seen, this quantity is sensitive to the details of the shape function. One 
can see this for oneself in Fig. |l^ for A = 0.3 GeV and Fig. [l7] for A = 0.5 
GeV. The contribution from the tree level term is presented in Figs. n3 and O 
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for a broader range of the hadronic mass. Similarly, the one- loop corrections 
themselves have been shown in Figs. [l4|, [l5[ 

5.6 Summary 

We calculate corrections of perturbative and nonperturbative nature to the 
hadronic mass spectrum in the semileptonic B decays. One-loop results are 
convoluted with the shape functions to give predictions of the \V u b\ matrix el- 
ement with an uncertainty of 10%. This error estimate relies on a the current 
knowledge of parameters characterising the B mesons. We state the functional 
dependence of the errors in terms of the involved parameters so that an esti- 
mate can be made in case the precision of the latter is improved on. We also 
suggest an observable to be employed in determinig these parameters from the 
semileptonic decays and present the prediction for it through one-loop level 
supplemented by leading twist corrections. 
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ntegrated hadronic mass spectrum — a s correction 
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Figure 9: Pcrturbativc correction to integrated hadronic mass spectrum. 
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Figure 10: Perturbative correction to integrated hadronic mass spectrum. 
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Convoluted ratios to order a =0.2 & 0.3 
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Figure 11: Ratios r for all the examined sets of parameters. 



6 Longitudinal polarisation of r lepton 
6.1 Introduction 

The longitudinal polarisation of the r lepton is a non-trivial quantity due to 
the large mass of the tauon. This is contrary to the electron or even muon, 
whose masses are so small that the V — A structure of the interaction causes 
them to be practically completely left-handed upon production in B decays. The 
polarisation itself is an interesting quantity as it can provide means of measuring 
the quark masses, especially the difference between b and c quark masses. Unlike 
the rate itself, the polarisation is regular througout the range of the charged 
lepton energy even after inclusion of one-loop corrections. Furthermore, the 
dependence on the poorly known CKM matrix elements is eliminated. The 
need for calculating the correction to this quantity is apparent given that the 
decay rate receives as much as 20% contribution from the one-loop diagrams. 
It is therefore all the more remarkable that the situation is in the end different 
for the polarisation. It turns out, as the calculation below shows, that the 
polarisation receives only a tiny correction. In view of the poor knowledge of 
the strong coupling constant, this kind of cancellation is also interesting when 
one considers extracting quark masses from this quantity. The first indication 
of the smallness of the radiative correction to the r longitudinal polarisation 
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Hadronic mass distribution — tree leve 
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Figure 12: Hadronic mass distributions at tree level. 
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Figure 13: Hadronic mass distributions at tree level. 
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Hadroric mass distribution — a s correction 
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Figure 14: 0(a s ) correction to hadronic mass distribution. 



Hadroric mass distribution — a, correction 



6 
4 
2 




A= 





3 GeV, X,= 


-0.1 GeV 2 


A= 





3 GeV, A,= 


-0.2 GeV 2 


A= 





3 GeV, X,= 


-0.25 GeV 


A = 





5 GeV, A,= 


-0.1 GeV 2 


A = 





5 GeV, A,= 


-0.4 GeV 2 


A= 





b GeV, A,= 


-0.6 GeV 2 



2.5 5 7.5 10 12.5 15 17.5 20 
hadronic mass squared in GeV 2 



22.5 25 



Figure 15: 0{a s ) correction to hadronic mass distribution. 
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Figure 16: Hadronic mass distributions to order a s for A = 0.3 GeV. 
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Figure 17: Hadronic mass distributions to order a s for A = 0.5 GeV. 
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was found in [[55] where this quantity was calculated in the rest frame of the 
pair of leptons. 

In order to find the polarisation, one may proceed as follows. First we get 
the result for the unpolarized distribution, as found in |3^, |37j. Then we are 
left with the task of finding the decay rate into either state of polarisation of 
the charged lepton. It seems natural, taking into account the structure of the 
weak interactions, that we calculate the rate of decay into a negatively polarized 
electron. Thus we present formulae for the double distribution in terms of the 
lepton energy and the invariant mass of the intermediate W boson. Combining 
these with the formulae for the unpolarized rate we readily find the polarisation. 
We also integrate the rate over the invariant mass of W to obtain the charged 
lepton energy spectrum. This allows us to evaluate the polarisation in terms of 
the tauon energy, which is shown in a plot. 

When giving the results for the polarisation we also use the nonperturbative 
corrections of order found in Q for the polarised lepton. That paper 

also gives the unpolarised distributions in the mass expansion, which were first 
calculated in |4|] and p3| . The HQET corrections appear in the form of extra 
terms proportional to the phenomenological constants Ai and A2. While the 
second of them can be related to the mass spliting between vector and pseu- 
doscalar B mesons due to chromomagnetic interaction, and is therefore fairly 
well established, the parameter Ai refers to the kinetic energy of the Fermi mo- 
tion of the b quark within the meson and only rough estimates can be made of it. 
It is thus convenient to express any corrections in terms of factors multiplying 
these parameters, which we do. This part of the paper starts with a description 
of the evaluation of the polarisation in Sec. 6.2, which is followed by the results 
given in an analytical form in Sec. 6.3. The method of incorporating the QCD 



corrections has been discussed in Sec. 3.1. Later on the moments of lepton 



energy distribution are studied, Sec, 6.4, and we summarize the results. 



6.2 Evaluation 

Let us start with the tree level matrix element for the unpolarized decay rate: 

MZ(r) = q ■ tQ ■ v. (144) 

With this kind of linear dependence on the four-momentum r, it is worth noting 
that the matrix element with the r polarisation taken into account is, 

K% = \mZ(K = r ma) - \{q ■ K)(Q ■ „), (145) 

where m stands for the lepton's mass and we have introduced the four-vector 
K 

K = t - ms. (146) 

This relation is rendered useful by expressing the polarisation four-vector s in 
the following way: 

s = At + BQ (147) 
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This is correct due to the fact that now only the temporal component of Q does 
not vanish, whereas the spatial parts of s and r are parallel. The coefficients A, B 
appearing in the formula above can be evaluated using the conditions defining 
the polarisation four-vector s: 

s 2 = -1 (148) 
s-t = 0. (149) 

Upon this one arrives at the following expressions: 

A ± = ±^-^—, (150) 

B ± = T JV!L. (151) 

r+ - r_ 

where the superscripts at A, B denote the polarisation of the lepton. 

Applying now the representation (147) of the polarisation s we readily obtain 
the following useful formula for the matrix element with the lepton polarized: 

= ^—^—M^ 3 (t) ± —H— M ^(Q). (152) 



The first term on the right hand side of (152) can be calculated immediately 
once we know the result for the unpolarized case. Thus the problem reduces to 
performing this calculation again with the only difference amounting to replacing 
the four-momentum r of the lepton with that of the decaying quark, Q. The 
Born-approximated r energy distribution can be written explicitly as 



UT f±(x), (153) 



dx 
where 

r o = ^|lW| 2 . (154) 
The Born level function fo(x) reads, for the unpolarized case, 

f (x) = ^( 2 t 3 {C[x 2 - 3x(l + r?) + 8 V ] + (3x - 6t7)(2 - x)} , (155) 

while the polarized cases are obtained using the function A/ : 

A/o(:r) = -^t 2 ( 2 {((3 - x - v )+3(x- 2)}, (156) 



in the following way: 



ti(x) = ^fo(x)±Af (x). (157) 



In the formulae above, 

r 3 = V* 2 - 47?, C = 1 - z P —— ■ (158) 
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At the tree level, one can express the polarisation integrated over the energy of 
the charged lepton as well: 



P = l-l/18{-24r; 2 + (4 l - 87 7 3 / 2 )(3 + r;-3p)+ 

I2 v (x 2 m ) - 3a4/2 + 3(1 - r,f{p - p 3 / s A ) + 

3(-l + n)p(l - p/s 2 )[3{l - v ) 2 + p(3 + 5r))\ + 3T 3 S/2 + 

3{x m - 2y/rj){-12r) - At] 2 + Ur/p - 3p 2 + 3 W 2 + p 3 ) - 

12 VP 3 Hs 2 /p) - 18(2r; 2 - p 2 - r, 2 p 2 ) In [2^j/(x m + T 3 )] + 

18(1 - 77 > 2 In (2,s 2 V^/[(l - rj){\ - v -T 3 )-p- rjp}) } 

/ {T 3 S/12 + {2 V 2 -p 2 - ifp 2 ) \n[{x m + T 3 )/(2V5j)] 

-(1 - tjV ln{Vw>/ [(1 - - V + T 3 ) - p - V p}}} , (159) 

where 

s=l~VV, T 3 = y/x 2 m - 47,, (160) 

S = l- 7[(1 + 77)(t7 + p 2 ) + p{\ + n 2 )] +f] 3 +p 3 + 12 VP . (161) 

Noting that the decomposition of the polarisation four-vector affects the lep- 
tonic tensor only, one sees immediately that it can also be applied to the QCD 
corrections. 



6.3 Analytical results 



The following formula gives the differential rate of the decay b — » tDX, X 
standing for a c quark or a pair of c and a gluon, once the lepton is taken to be 
negatively polarized: 



dT- 
dx dt 



i2r n 



Fo(x,t) 



Fr A {x,t) 



3?r 1,A 



i2r f^- B (x,t) 



for (x, t) in A 
for (x, t) in B 



(162) 



F^~ differs according to which region (A or B) it belongs to. Region A is available 
for the 3- and 4-body decay, while Region B with a gluon only. The following 
formulae are given for the negative polarisation of the lepton, that is, we take 



A- = - 
B~ = - 



1 



2VV 



(163) 
(164) 



The factor Tq is defined in Eq.(154), while 



F- (x, t) = (1 9 X + t) [t + (x -t- V )- V (l+p- t)] (165) 
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and 



F^ A (X, t) = F -$ + D n*n + D 
n=l 
5 

F- B {x, t) = F -* + J2 D%*n + D, 



.1 

6 > 



(166) 



(167) 



71=1 



The factor 12 in the formula (|162| ) is introduced to meet the wid ely used |55|, 
p5[ |59| convention for Fo(x) and To. The symbols present in (16£) are defined 
as follows: 



2Po , , 1 - r H 

*o = -Ni 

P3 P+ 



-£i 2 (l- 



1-^ 



ii 2 (l 1—) 

P+ 

l-t/r+ 
Li 2 (l ' 

P- P- 

+Li 2 ( w- ) - Li 2 ( w+ ) + AY p In y/p\ 
- p) -Ah\z maxi 



+4(l-_-y p )ln(z 

P3 



$i = Li 2 (w-) + Li 2 (w + ) - Li 2 (r + ) - Li 2 (t/r + ) 

$2 = — , 



P3 



1 



*o = 4(^Kp-l)ln( 

P3 



2 P°IT ■ ft 

[Li 2 (l - 

P3 

-Li 2 (l- 
+Li 2 (l- 
-Li 2 (l - 



*3 = -ln^/p, 
$ 4 = -ln(l-r + ), 

$ B = 5Mi-t/T+) 1 

) +41n( ) 

1 - T +-) + Li 2 (l- 1 -^t±) 



P- 

1-T+ . 
1-T- . 
1-T-, 



P- 



£»2(1 



1 - t/r- 



P- P- 
*i = Li 2 (r+) + Li 2 (t/r + ) - Li 2 (r_) - Li 2 (t/r_ 

* 2 = - ln(l - t_ 
* 3 = iln(l-t/r_ 
*4 = - ln(l - r+ 



(168) 
(169) 
(170) 

(171) 
(172) 
(173) 



(174) 
(175) 
(176) 

(177) 

(178) 
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* 5 = iln(l-t/r+). 



(179) 



We introduce C1...C5 to simplify the formulae for and D 



Ci = 



1 



[p(— rfx + Arjt — rjT + — 477 + xt + — tr+) — 2p 2 rj + rj{xt + x 



+At - 2t 2 + 6 + 27?) +T+{--qx + rft + 5r? - 2xt + x + t + t 2 )] ; (180) 
C 2 = p(r] - t - 2t+) + (xt 2 -2t 2 T+)/(2r 1 ) +j](-x -2t + 2T+ + 30)/2 



+p 2 (-9rj + 7+) + (2xt 2 T+ - x 2 t 2 /2 - t 2 rl)/ V + 77 2 (-t+ + 10) 
+xt(x - t - l)r](4xt - xt + -2x- x 2 /2 + tr+ + At 

-At 2 + llr+ - t 2 - 12) + +r+(-6xt - x + 2x 2 - 3t + 5t 2 - 5)] ; (182) 

C 4 = . 1 \p(2rjXT + /t - 8-qx + 8r]T + /t + 18r)t - II777+ - 4?7 2 T + /t 2 
y/x 2 - Ar\ 

-xt + 5t + x - 7tr + ) + p 2 (2r]T + /t - 9 77 - f] 2 T + /t 2 ) + xt{\ +t-x) 

+ {-2xt 2 T + + x 2 t 2 /2 + t 2 Tl)/r] + r, 2 {hT+/t 2 - 6r+/t - 16/t + 6) 

+r](2xT + /t + 2xt - 3xt+ -Ax + x 2 /2 - 10r+/t + 3tr + + 12t 

-At 2 + 3t+ + rl + 4) + r+i-Axt - 3x + 2x 2 + lit + 2t 2 )] ; (183) 

C 5 = [ P /(x-t-r ] /t)}(-p7 1 2 T + /t 2 + f ] T + /t-f ] -f ] 2 T + /t 2 +7 1 2 /t + pf ] T + /t 

-m + PV 2 /t)/2 + [p/(l -x + 7j)](-7?i - 7/T+ + 77 2 + tr+)/2 
+p{3r]T+/t + 977 - 9t - 3t+)/2 + (xt 2 /A - t 2 T + )/rj + r)(3x - 10r+/i 
-12* - 2t+ - 24 + 2t7)/4 + (3t + 5)t+/2 - xi + 6t + 5t 2 /2; (184) 
Df = d; (185) 

= — 7 ^^={pr](3A-6xT + /t 2 -Ax/t 2 + 3x/t + 5xt-2xr + + 2x 
A^Jx 2 - Ar] 

+3x 2 /t 2 + x 2 + 8r + /t 2 + 18r+/i - 2/t + 2tr + + 38t - lAt 2 + 20r+) 
+p7/ 2 (16 - 2x/t 2 - hr+ft 2 - 2r + /t + 16/t - 7+) + pr + (-Axt - Ax + 6t 
+7t 2 + 11) + p 2 [t?(-10 + 6xr + /t 2 + 6x/t 2 + Axr + /t - 3x - 3x 2 /t 2 
-2x 2 /t - 12r + /t 2 - 18r+/i + 6/i + 18f - 6r+) + 77 2 (a;/t 2 + r + /t 2 
-r+/t - 2/t) + t + (2x - 5t - 7) + +p?7(-10 - 2xr + /t 2 - Ax/t 2 - x/t 
+x 2 /t 2 + 8r+/t 2 + 6r+/i - 6/i) + pv 2 r+/t 2 + pr+ + p 2 ri{x/t 2 - 2r+/t 2 
+2/t)] + T + (-Axt + 2xt 2 + 2x + 7t + t 2 - 3t 3 - 5) + 7?(-14 + 2xr + /t 2 
+x/t 2 - Axr + /t - 2x/t + Axt - 2xt 2 + 2xt + -x- x 2 /t 2 + 2x 2 /t - x 2 
-2r + /t 2 - 6r + /t - 10tr + + 32t - 22t 2 + At 3 + I87+) + t/ 2 (28 + x/t 2 
-2x/t +x + 3r + /t 2 - 5r+/i - U/t + tr + - Ut + r+)} ; (186) 



C 3 = 




(181) 
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Df = {p [q(-6 - 4xr + /t 2 - 3x/t 2 - 2xr + /t + x/t - 13x + 2x 2 /t 2 

y?X 2 — At] 

+x 2 /t + 6r + /t 2 + 12r+/i - 2/t + 26t - 10r+) + r] 2 (-x/t 2 - 2r+/t 2 
+2/t) +t + (6x + 6 - 12i)] + p 3 r)(-x/t 2 + 2r + /t 2 - 2/t) 
+p 2 [»?(-10 + 2xT + /t 2 + 3x/t 2 - a; 2 /* 2 - 6r+/t 2 - 4r+/t + 4/t) 
-ri 2 T + /t 2 - t+] + ?7 2 (14 + x/i 2 - x/t + 3r + /t 2 - 2r + /t - 14/i - r+) 
+7?(-14 + 2xr + /t 2 + x/t 2 - 2xr + /t - x/t + 2xt -2x- x 2 /t 2 + x 2 /t 
-2r + /t 2 - 8r+/t + 18t - 4t 2 + 10r+) + r+(-2.xt + 2x + 2i + 3t 2 - 5)} ; 

(187) 

Dt = -C 3 -C 2 ; (188) 
Dt = -C 4 + C 2 ; (189) 

Dt = -\c 5 + - T= !=={[p/(l-x + V )}[ V (tT + +t + T + )-r 1 2 (l+t + T + ) 
* 4a/ x z — 4?7 

+ V 3 - tr+)] + [p/(x - t - v /t)](l + p) [v(t - r+) + r/ 2 (-l + r+/t 2 + 

T+/t - l/t) +Tf{l/t 2 - r+/t 3 )] + pry(43 - arr + /i - Ax/t - 5x + 2x 2 ft 

+9r + /t + 13t + 5r+) + pr) 2 (-5 - r + /t 2 + 2r + /t + l/t) + T + (9xt + 5x 

+At - 6t 2 ) + 12xt + hxt 2 - 2x 2 t + p(-9xt - 3xr+ + 5tr+) + 

p 2 T](-lO + 2x/t - 3r+/t) + p 2 rf{-T + /t 2 + l/t) + {-6xt 2 T + + x 2 t 2 

+2t 2 r 2 )/(2ry) + r/(-35 - xr+/t + 2x/t - 2xt + Axt+ + 26x 

-2x 2 /t - 3x 2 /2 - Ar+/t - 2tr+ - Mt - t 2 - 22t+ - 3r 2 ) 

+7] 2 (-2 + 2x/t-x + 6T+/t + 2t)} ; (190) 



D? 


= c x 


(191) 


Di 


= C2 — C3 


(192) 


D§ 


= — C2 — C4 


(193) 


Df 


= C2 + C3 


(194) 


Di 


= — C2 + C4 


(195) 


Dg 


= C 5 


(196) 






(197) 



One can perform the limit p — > 0, which corresponds to the decay of the bottom 
quark to an up quark and leptons. The formulae, which are much simpler in 
this case, are presented in the same manner as the full results: 



d£_ 

dxdt 




F^{x,t)-^-Fr A {x,t) 



o^rF- B (x,t) 



for(x, t)inA 
for(x, t)inB 



with 



F~(x,t) 



(l-x-t) 

T+ -T_ 



[r+(x - t - V ) - V (l - t)} 



(198) 



(199) 
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(I- 



F- A (X, t) = F $o + Y, V n®n + V, 

71=1 

5 

F- B {X, i) = F f + ^ V n^n + Vf. 



(200) 
(201) 



n=l 



*o = 2 



Li 2 (^-4) + Li 2 { l ! T + 1 ) + Li 2 (t) 



1 



l-t ' 1/t- 1 

+ ln 2 (l - t+) 2 ln 2 (l - t) + ln 2 (l - t/r+) 
-21n(l - t) lnz maXl 

ir 2 

*i = 777 + Li 2(t) - Li 2 (T+) - Li 2 (t/T+), 



12 



2o3 



21n(l-t) 
l-t ' 

$4 = $ 4 , 
$5 = $5, 



(202) 

(203) 

(204) 
(205) 



*n = 2 



L *2(^— r) + £*2(-r7T — 7-) - Lz 2(-r— r) 



1 -i 7 ' " v l/t - 1 

-1/7+ "I 



-Li 2 (- 



1/t — 1 



+ ln(l - t) ln( 



l-t 

z max > 



-ln(i-^)ln[(l-r + )(l-r_)] 

■In(^^)In[(l-t/r + )(l-t/r-)], 
*„ = *„, (n=1...5). 



(206) 
(207) 



a/x 2 — 4 77 
T + (x-2xt + t + t 2 )] 



[r?(6 + xt - xt+ + x + tr++At- 2t 2 + 5r + ) + 2rf 



(208) 



{xt 2 - 2t 2 T + )/(2i]) + 7/(30 - x - 2t + 2r+)/2 - xt - 2xt+ + 2tr+ 

+t + t 2 + 6r+; (209) 

[(4xt 2 r+ - x 2 t 2 - 2t 2 r 2 + )/(2 V ) + + V 2 (1Q - r+) 

y/X Z — 4?7 

T](-24 + 8xt - 2xt + - Ax - x 2 + 2tr + + 8t - 8t 2 + 22t+ - 2t 2 )/2 



50 



-6xtr + -xt- xt 2 - xt+ + x 2 t + 2x 2 t + - 3tr + + 5t 2 r + - 5t+] ;(210) 

C£ = 1 [(~4xt 2 T + +x 2 t 2 + 2t 2 Tl)/(2f])+r](4 + 2xT + /t + 2xt 

x 2 — 4?y 

-3xt+ -Ax + x 2 /2 - lOr+ft + Mt+ + I2t - At 2 + 3r+ + r|) 

+?7 2 (6 + 5t+/< 2 - 6r + /t - 16/*) - 4artT+ + art + xt 2 - 3iet + 

-£ 2 * + 2x 2 t + + 11*t + + 2* 2 t + ] ; (211) 

= (x* 2 -4* 2 t+)/(4?7)+7 ? (-24 + 3:e-10t + /*-12*-2t + )/4 

+r/ 2 /2 - art + 3*r+/2 + 6* + 5t 2 /2 + 5r+/2; (212) 

= d (213) 

D^3 = 1 [r/(-14 + 2xr + /t 2 +a;/t 2 -4xr + /t-2a:/f + 4xf-2xf 2 

4^x2-477 

+2iet + - x - x 2 /t 2 + 2x 2 /t - x 2 - 2r + /t 2 - 6r+/* - 10*t+ + 32t 
-22t 2 + 4t 3 + 18t+) + r] 2 (28 + x/t 2 - 2x/t + x + 3t+/* 2 - 5r+/t 
-14/* + tr+ -14* + r+) + (-4artT + + 2x* 2 t + + 2xt + + 7tr + + t 2 r + 
-3* 3 r+ -5r+)] ; (214) 
Vf = -Ca-C 3 , (215) 
^ = C 2 -C 4 , (216) 

2^ = -k, + J_= [(-6x* 2 t + + x 2 t 2 + 2t 2 r 2 )/(277) 
2 4ya; 2 - 

+?7(-35 - :et + /* + 2x/t - 2xt + 4.tt+ + 26a; - 2x 2 /t - 3x 2 /2 - 4t+/* 
-2tr + - 34* - t 2 - 22t+ - 3t 2 ) + r\ 2 {-2 + 2x/t -x + 6r+/* + 2t) 
+ {9xtr + + Ylxt + 5xt 2 + 5xt + - 2x 2 t + Atr + - 6* 2 r+)] ; (217) 

On integration over t, one obtains tau energy distributions according to the 
formula 

121 o dx Sir 

where if we drop the superscript '— ' we obtain the corresponding formula for 
the unpolarized case. 

The results are presented in Fig. [l8|, where the polarisation is plotted versus 
the t lepton energy. Both the Born and first order approximation are showed. 
In order to make the correction more explicit, we also present another diagram, 
where the function R(x) is drawn, defined as following: 

i-jw-p-ftwm + fij*,)] (2i 9) 

which gives the expression for R(x): 

*m = W\-tt\' (220) 

fo(x) f (x) 
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Polarization of the r lepton 
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Figure 18: Polarisation of r lepton in the Born approximation (dashed line) and 
including the first order QCD correction (solid line) as functions of the scaled 
t energy x. The mass of the b quark taken at 4.75 GeV, c quark 1.35 GeV and 
the coupling constant a s = 0.2 

thus its meaning is how the radiative corrections differ with respect to the state 
of polarisation they act on. In ( |219|) Pq denotes the zeroth order approximation 
to the polarisation. The function /£" (x)/fg(x) has been presented in Fig. [l9], 
while Fig. ^ shows the function R(x), so that one can immediately see how 
small the correction is in comparison to the correction to the energy distribution. 

Integrating over the charged lepton energy one can obtain its total polarisa- 
tion as well as the corrections to which it is subject. If we take m& — 4.75 GeV 
as the central value for the decaying quark mass and mb = 4.4 and nib = 5.2 for 
the limits, we arrive at the following (the mass difference mj — m c — 3.4 GeV 
everywhere) : 

l-P-d-A) + + + (221) 
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QCD correction for polarized t 
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x = 2E r / n b 



Figure 19: The ratio f l (x)/ f (x) representing the radiative correction for the 
negatively polarized state as dependent on the scaled r lepton energy x. 



with 



Pa - -0.7388 + Q;°i°g 

R s = — 0.016 + q;oi7 

k-1 - 421+ 027 

E>2 n OQ — 0.22 



6.4 Moments of r energy distribution 

The moments of r energy distribution, which are useful sources of informa- 
tion on the physical parameters regarding the discussed decay can be evaluated 
according to the formula: 



Mi 



(222) 



± 



M 

rt = (223) 
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Figure 20: The QCD-correction function R(x) for the pole mass values of the 
b quark taken to be 4.4 GeV (dashed), 4.75 GeV (solid) and 5.2 GeV (dash- 
dotted) as dependent on the scaled r lepton energy x. 



where E m i n and E ma x are the lower and upper limits for r energy and M n 
include both perturbative and nonperturbative QCD corrections to r energy 
spectrum. The superscripts denote the polarisation states. Since one obviously 
has 

M n = M+ + M~ (224) 

where M n stands for the unpolarized momenta, we only give the values of the 
momenta for the negative polarisation case. The unpolarized distributions were 
given in p6| , |37| . The nonperturbative corrections to the charged lepton spec- 
trum from semileptonic B decay have been derived in the HQET framework up 
to order of 1/ m\ ||^, |i[ [f5| . The corrected heavy lepton energy spectrum 
can be written in the following way: 



i dr 




(225) 



where Ai and A2 are the HQET parameters corresponding to the b quark kinetic 
energy and the energy of interaction of the b quark magnetic moment with 
the chromomagnetic field produced by the light quark in the meson B. The 
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functions /np' 2 can be easily extracted from the formula (2.11) in |f45f.The 
formula ( |225| ) looks identically if one considers definite polarisation state of the 
final r lepton. The appropriate calculation within the HQET scheme has also 
been performed [ft5f , see formula (2.12) therein. 
Following ref. |59(| we expand the ratios r n : 



V = T 



+ + (226, 



where is the lowest approximation of r„, 

Each of the 6n is expressed by integrals of the corresponding correction func- 
tion fw(x) and the tree level term fo{x) 

£%- p ^-w* n *° s^rf^ (228) 



where the index i denotes any of the three kinds of corrections discussed above. 
The coefficients Sn depend only on the two ratios of the charged lepton and 
the c quark to the mass m b . Following ref. [fj6[ |3?j] we employ the functional 
dependence of the form 

5^-(m b ,m c ,m T ) = ^) . (229) 

The quark masses are not known precisely so we have calculated the coefficients 
in a reasonable range of the parameters, that is, 4.4 GeV< rrib < 5.2 GeV and 
0.25 < m c /mb < 0.35 and then fitted to them functions of the following form: 

S(p,q) = a + b(p - p ) + c(q - q Q ) + d(p - p ) 2 

= +e(p-Po)(q-qo) + f(q-qo)\ (230) 

where p = m b /m Tl p = 4.75 GeV/ 1.777 GeV= 2.6730, q = m c /m b ,q = 0.28 
and the polynomial coefficients can be fitted for each of the Sn separately 
with a relative error of less than 2%. Our choice of the central values reflects 
the realistic masses of quarks: m b = 4.75 GeV and m c = 1.35 GeV, for which 

5 {l) - a W 

To bring out the difference in the extent to which the corrections affect the 
two different polarisation states it is useful to compare these coefficients with 
the ones obtained with the polarisation summed over. This can be done along 



the lines suggested by the treatment of the polarisation itself, see Eq.(219). The 
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corresponding expansion takes the form 



-(0) 



„(0) 



2a s 
3tt 



(231) 

With these, one can readily find the actual relative correction of each kind, 
assuming reasonable values of a s , Ai andA 2 . Here we take a s = 0.2, 0.15 GeV 2 < 
-Ai < 0.60 GeV 2 , A 2 = 0.12 GeV 2 , keep the b quark mass fixed at 4.75 GeV 
and the mass of the c quark equal to 1.35 GeV. The corrections then read for 
n = l(n = 5) 





Correction to 


Type 


- / o - 


/ (o) 


perturbative 


-0.00084(-0.0048) 


-0.0009(-0.0052) 


kinetic energy 


0.008 ±0.005(0.06 ±0.04) 


0.008 ± 0.005(0.06 ±0.04) 


chromomagnetic 


-0.0097(-0.053) 


-0.0092(-0.0511) 



6.5 Summary 

We have calculated the one-loop QCD corrections to the longitudinal polar- 
isation of the r lepton in semileptonic B decays. This calculation has been 
published in jl9]]. As was expected from earlier results (3f|, the effect of the 
correction on the polarisation is negligible, which makes it a promising quantity 
for the determination of the quark masses. Discussing the moments of the lep- 
ton energy distribution, we have also included the nonperturbative corrections. 
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7 Polarisation of r lepton with respect to W 



7.1 Introduction 

The method used in the calculation of the longitudinal polarisation can easily 
be modified to give other polarisations. This fact matters insomuch that ex- 
perimentally it is the polarisation along the intermediating W boson direction 
that is easier to measure [Q , see also ||7J . The reason is that the direction of r 
lepton can be determined at B factories with rather poor accuracy. On the other 
hand the direction of W is opposite to the direction of hadrons in semileptonic 
B decays. The latter can be well measured at least for the exclusive B — > Dtv t 
and B —* D*tv t channels which probably contribute the dominant contribution 
to the inclusive decay rate. 

In the calculation below, we find the tree level term and the 0(1/ m^) HQET 
corrections to the doubly differential rate with the final r lepton polarized along 
the momentum of the intermediate W boson. The one-loop correction is pre- 
sented in a plot as it has been evaluated numerically. Also we give the tree level 
charged lepton distribution in an analytic form. 

The HQET mass power corrections to the lepton spectra in semileptonic B 
decays have been applied to the calculation of the polarisation of the lepton. The 
mass of the lepton had of course to be taken into account, which was done for 
unpolarised distributions in |ff2| , [f3[ . The kinematical complications arising due 
to the massive lepton are also discussed therein. The longitudinal polarisation 



was calculated in 45 1. Some further problems emerged in the calculation of the 
corrections to the transverse polarisations found in p8| . The case we consider 
leads to apparent divergences, but they cancel so that no special regularisation 
has to be employed. 



After a brief note on the kinematics in Sec. 7.2 we discuss the way the po 



larisation is evaluated, Sec. 7.3. An analytic formula for the double differential 



distribution is presented including the mass power corrections in Sec. 7.4. Due 



to its complicated form, we only show numerical form of the lepton energy dis- 



tribution in Sec. 7.5. This section contains also the one-loop perturbative cor- 
rections. All of these corrections are presented in a plot. Finally we summarize 
the results. 

7.2 Kinematics 

Recall that the phase space is defined by the ranges of the kinematical variables: 

2y/T] <x< l + rj-p = x m , (232) 

tmm = T- f 1 - 1 _ P T ^J - t - T + f 1 - Y^V) = tmaX ' ^ 233 '' 

The limits above are obtained within the parton model approximation. They 
change if we allow for Fermi motion, which we must in order to be able to 
discuss the HQET corrections to the decay widths |1[ |§ f|| g g Also, 
contrary to the parton model case, the energy of neutrino can vary within limits 
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which depend in a non-trivial manner on the values of the variables x, t. The 
details of the subject were discussed in E^], so we will only state here that the 
integrations involving delta functions and their derivatives have the effect of 
confining the range of the variables x, t to that of the parton model. This is 
easily seen from the form of the corrections, see Eq. (p2|) . 

7.3 Polarisation evaluation 

The polarisation is found by evaluating the unpolarized decay width and any of 
the two corresponding to a definite polarisation, according to the definition, 

r+ - r~ r~ 
P = F+TT^ = 1 - 2 - ' (234) 

where T = T + + T~ . In this calculation we follow the method employed when 
dealing with the longitudinal polarisation of the charged lepton. Thus in the 
rest frame of the decaying quark, one can decompose: 

s = AQ + BW . (235) 

The coefficients A, B can be evaluated using the conditions defining the polari- 
sation four- vector s, which reduce to the following when the parton model value 
of the neutrino energy is assumed: 

At = T ; t + V , (236) 

\/t(x — X-)(x+ — x) 

X 

y/t(x — x-)(x+ — x) 
where the superscripts at A, B denote the polarisation of the lepton, while 

x± = (l + r)/t)w± . (238) 



B± = ± j= 1= == , (237) 



However, the coefficients in the decomposition (147) need to be re-evaluated, 
taking into account the Fermi motion and working in the rest frame of the 
decaying meson. We will presently construct the four-vector s representing a 
charged lepton polarized along the direction of the W boson. The defining 
properties of s are 

s 2 = -1 (239) 

and 

s ■ t = , (240) 

complemented by the relation s \\ W. Since we are working in the rest frame of 
the decaying meson, the four-vector s can be decomposed as 

s = Av + BW , (241) 

where v and W stand for the four-velocity of the B meson and the four- 
momentum of the intermediate W boson, respectively. While this form au- 



tomatically satisfies s \\ W, the other two relations (239) and (|240|) have to be 
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imposed, hence yielding the expressions for the coefficients appearing in (241). 
With v = (1, 0, 0, 0), one readily identifies: 

vt = x/2, vv = x l/ /2, v 2 = -1 . (242) 

These combined with the other dot products lead to the following formulae: 

A ± = T{t + 71) , (243) 

y/x(t + rf)(x + x u ) - xH - (t + i]) 2 

B± = ±X = . (244) 

y/x(t + rf)(x + Xu) - X 2 t ~(t + T]) 2 

The evaluation of the HQET corrections involves differentiation over the neu- 
trino energy, once or twice. The denominator in the above expressions is easily 
seen to vanish at the point where the W boson is at rest, ft is known that the 
kinematics of the process, together with the delta functions and their deriva- 
tives, finally reduces to integration over the partonic phase space. Then there 
is one such point where the denominator vanishes, 

x = l-^+_?L_ ( (245) 

* = (1-Vp) 2 • (246) 

One might thus raise the question of analyticity of the expressions obtained 
in this way. However, the divergences cancel and moreover the resulting dis- 
tribution is continuous if we ignore the endpoint behaviour. That this is so 
indeed, may be verified by changing the variables from t to the square of the 
three-momentum of the W boson. Then the singularity makes its presence 
only on integration over w| rather than affecting the analytical structure of the 
distributions. It turns out that using this variable one obtains an analytic ex- 
pression. This is made clear once one notices that the only terms that occur in 
the course of the calculation are the dot products of the four-vector s and the 
other four- vectors. Writing them out explicitly, 

cos 9 

s+-v= -, a 2 , , (247) 

\Jx z — T3 cos^ 6 

^{x + x^r^e-^ 

2Vx 2 - t 3 cos 2 e 

with 

cosg = W 3-^-(^-^)/4 (24Q) 
W3T3 

where the subscript denotes the polarisation direction, we easily verify that 
the triple differential distribution is analytic in the neutrino energy. Lastly, 
let us note that another change of variable can be useful for evaluating the 
distribution. Namely, using the cosine of the angle subtended by the tau lepton 
and the neutrino |4q] eliminates the singular terms from the double differential 
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distribution. We have checked numerically that the resulting distribution is the 



same. Applying now the representation (235) of the polarisation four- vector s 
we obtain the following formula for the matrix element with the lepton polarized: 



Im uu (t) ± ^ 

2 yjx(t + r)){x„ + x) - X 2 t ~(t + l]f 

-(t + r 1 )M un (Q)} . 



[xM un (W) 



(250) 



The above expression is valid for the HQET corrections, too. The first term 



on the right hand side of (250) can be calculated immediately once we know 
the result for the unpolarized case. Then the other terms require the formal 
replacement of the four-momenta r — > W and t — ► Q in the argument. 



7.4 Double differential distribution 

The polarized distribution can be written in the form, 



To dx dt 2 

where 



' ' /] ]-F un ± ( F + F+ - F_ ) , (2.-)l) 



r ° = W IW|2 - (252) 



The first term on the right hand side of Eq.(251) stands for the unpolarized 
distribution, given e.g. in [j42f , Eq.(30). Here we will only give the new other 
term: 1 

F — y/rj W [6/r + K b W (/ 2 + f 3 W 2 + |/ 4 W 4 ) + G b (/ 5 + f 6 W 2 )] (253) 
where 

/i = -xt(l + p - rj) + 2xt 2 - x(l + pi] - 2p + p 2 + T]) - x 2 t 

+x 2 {l -p)+ t(l + 2pri - p 2 ) + t 2 (2p - rj) - t 3 - p 2 V + r] , (254) 

f 2 = -8xt(l - p + rj)- IQxt 2 + 8x(l + pj] - 2p + p 2 - rj) + 6x 2 t 

- 6x 2 (l - p) + 2i(l - 6pr] -Ap + 3p 2 + 12?/) + 2t 2 (8 -6p + 3rj) 
+ 6t 3 - 8p?7 + 6p 2 r] + 2r] + 8r] 2 , (255) 



f 3 = xt(- 10 prj + 6 prf + Up 2 rj + 9p 2 if - 6p 3 r] + 2ij + 9rf - 4rj 3 ) 
+ xt 2 (l - 9/V - 5p + I8p 2 -q + 7p 2 - 3p 3 + 6?7 - 13r/ 2 ) 
+ xt 3 (l - 18/9?7 - 2p + 9p 2 - 1477 + 3t? 2 ) - xt 4 (5 + 9p- 677) 
+ 3xt 5 + x{~bprf + Apr] 3 + 7p 2 rf - 3p 3 i] 2 + rj 2 + Arf) 
+ x 2 t(l + 14/97/ + 2 P T] 2 - 3p - p 2 r[ + 3p 2 - p 3 - 13tj + 8r/ 2 ) 
+ x 2 t 2 {-6 + 5p77 + 6/9 + 1577 - ?7 2 ) + x 2 t 3 (7 + 3p- 3??) - 2x 2 t A 

X A FORTRAN code for this formula is available from piotr@charm.phys.us.edu.pl 
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+ x 2 r)(-3p + 8/977 + 3p 2 - p 2 r] - p 3 - 7t? + 1) + x 3 t(l - p) 2 - x 3 t 2 rj 
-x 3 t 3 +x 3 r](l - p) 2 , (256) 
fi = x 2 t(-6pr] 2 - 2p 2 rf + 6p 3 r/ 2 + Ap 3 rf - 3/?V + 3t? 2 - 2rf) 

+ x 2 t 2 (-6pr] - 2prf - 6/?V - Qp 2 rf + 6p 3 r/ + 12pV - 3/9 4 t? + 3?? - 6r/ 2 ) 

+ x 2 t 3 (l - Qprj 2 + Apr, 3 -2p- 6p 2 r/ - 18/9 2 t? 2 + \2p 3 r] + 2p 3 - p 4 - 6t? + 2t? 3 ) 

+ x 2 t 4 (-2 - 6/977 + \2prj 2 - 18p 2 rj - 2p 2 + Ap 3 + 6t? 2 - rj 3 ) 

+ x 2 t 5 (l2pri -2p- 6p 2 + 677- 3t? 2 ) + x 2 t 6 (2 + Ap - 3??) - x 2 t 7 

+ x 2 (-2pr] 3 + 2p 3 rf - p 4 r/ 3 + rj 3 ) + x 3 t(8pr/ + prf 2 + 2pr/ 3 - 12p 2 r/ + p 2 r/ 2 

+3p 2 rj i + 8p 3 r] - p 3 rj 2 - 2p 4 r/ - 2r] - rj 2 + 3r] 3 ) 
+ x 3 t 2 (-l - pr} + 6/9T? 2 - 3prf + Ap- p 2 rj + 9p 2 rf - 6p 2 + p 3 i] + Ap 3 

-p 4 + V + 9r/ 2 - 3t? 3 ) 
+ x 3 t 3 (l + 6pr] - 11/97? 2 - p + 9p 2 r] - p 2 + p 3 + 9rj - llr] 2 + T] 3 ) 
+ x 3 t 4 (3 - 13/37? + 2p + 3p 2 - 13/7 + Arf) + a; 3 7j 5 (-5 - bp + brf) + 2x 3 t 6 
+ x 3 (Apr/ 2 + prj 3 - 6p 2 r] 2 + p 2 r\ 3 + Ap 3 rj 2 - p 3 i] 3 - p 4 r/ 2 - rj 2 - rj 3 ) 
+ x 4 t(-6pi] + 2pr/ 2 + 6p 2 r/ + p 2 rf - 2p 3 rj + 2i] - 3t? 2 ) 
+ x 4 t 2 (l + Apr] + prf -3p + 2p 2 r] + 3p 2 - p 3 - 677 + 3rf) 
+ x 4 t 3 {-3 + 2p V + 2p + p 2 +6ri- rf) + x 4 t 4 (3 + p-2rj) 

- x 4 t 5 + x 4 {-3prf + 3p 2 r l 2 - p 3 rf + rf) , (257) 
f 5 = Axt(3 + bp- 5j?) - AOxt 2 + Ax(l + bpr/ -6p + bp 2 + rj) + lQx 2 t 

- 2x 2 (l - bp) - 6t(l + 10/97? - bp 2 ) - 2t 2 (A + 30p - 15t?) + 30t 3 

+ 30p 2 ri - 6t? + 8?? 2 , (258) 
f 6 = xt(-l8prj - 2pri 2 + 6p 2 r] - lbp 2 r/ 2 + 10p 3 r/ + 2r/ + 17r? 2 - Arf) 
+ xt 2 (l - 16/97? + Ibprf - 9/9 - 30p 2 r/ + 3p 2 + bp 3 + 10t? - t? 2 ) 
+ xt 3 (l + 30/97? - lOp - lbp 2 + 10t? - 5t? 2 ) + xt 4 (7 + lbp - 10r?) 

- bxt 5 + x(-9pr/ 2 + Apr/ 3 + 3p 2 r/ 2 + bp 3 r/ 2 + rf + 8r/ 3 ) 

+ x 2 t(-l + 30/97? - 10/97? 2 + 7p + bp 2 r/ - lip 2 + bp 3 - llr? - 4t? 2 ) 
+ x 2 t 2 (-2 - 25/jt? + 18/9 - 19r? + brj 2 ) - 15x 2 t 3 (l + p - rf) 
+ 10x 2 t 4 + x 2 (7pr] + 12/97? 2 - llp 2 r/ + bp 2 r/ 2 + bp 3 r/ - 7? - 9t? 2 ) 
+ x 3 t(l - 6/9 + bp 2 + 8t?) + x 3 t 2 (8 - 5t?) - bx 3 t 3 

+ x 3 (-6/9?? + 5/9 2 7? + 7/) , (259) 

F± = W± { [K h (h li± + h 2 ±W± 2 ) + G h h 3 , ± ] 5(z±) 

+K b h 4 ,±5'(z±)} , (260) 

where 

h h ± = -8xtr/ + 8xri 2 + Ax 2 t + 12x 2 r) + 2a; 3 t - 2x 3 r\ - 2x 4 - 16£t? - 16r? 2 
- 4cr± (6xt - 2xrj + 3x 2 t + x 2 7? - a; 3 - 8t-q - At 2 - Arj 2 ) 
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- 8(t 2 ± (3xt -xr]- bx 2 + At + Arj) + I6a^(3x + t + rj) , (261) 
h 2 ,± = 2a± x (8tr] 2 + 16t 2 r] + 8t 3 - Axtrj 2 + Axt 3 - 8x 2 t-q - 6x 2 t 2 - 2x 2 rr 2 

- x 3 t 2 + x 3 r} 2 + xH + x A rj) + A<7 2 ± x [-Atr] 2 - 8t 2 r] - At 3 
-\2xt-q - 8xt 2 - Axrf 2 + 2x 2 tr] - x 2 t 2 + 3a; V + 3x 3 t + 3x 3 rj) 
+ 8cr± x (-Atr) - 2t 2 - 2r\ 2 + Axtrj + xt 2 + ixrj 2 + 3x 2 t + 3x 2 rj) 

+ 16<4 x (2tr) + t 2 + -q 2 + xt + xt]) , (262) 
ha,± = -I6xtr) - 8xt + 8xt 2 + 8xr/ + 8xr] 2 - 8x 2 t + 8x 2 r] + 2At 2 -2Ar} 2 

+ 4cr± (-10xt - Uxr) - 5x 2 t + 5x 2 r/ - Ax 2 + 5x 3 + I2t - At 2 + 12?y 
+ Ar/ 2 ) + 16ct| (bxr] -2x + bx 2 - 9t - 9i]) 

+ 80a 3 ± (x + t + rj) , (263) 
/i 4 ± = Aa± {Axtrj - Axt 2 + 6x 2 t + 2x 2 r) + x 3 t - x 3 rj - x 4 - 8tr) - 8t 2 ) 
+ 8a± (8xt + Axr] + x 2 t - 3x 2 r) - 3x 3 + Atr\ + At 2 ) 

A 



16a% (xt + 3xrj + 3x 2 - 2t - 2rj) - 32ai (x + 1 + 77) , (264) 



and 



W± = - , (265) 

y/ X (t + 77)(20-± +X)~ X 2 t - (i + 7 ? ) 2 

W = 1 , (266) 

y/x(t + r))(x v0 +x)- xH - (t + rj) 2 

a± = (t - rj)/(2r±) , z± = l + t- p-x-2cr± . (267) 

The parameters Kt,, Gb, representing the kinetic energy and the chrome-magnetic 
energy, are defined according to 0] 

7.5 Lepton energy distribution 

As regards the HQET correction terms, we only give the energy distribution in 
the form of a diagram evaluated numerically. Beneath we also give the Born level 
approximation analytically. The analytic formulae for the polarized distribution 
can be simplified if we split the kinematical range of y into two parts, separated 
by the value of the charged lepton energy where the virtual W boson can stay 
at rest. This value is 

In the formulae below, the superscripts A, B refer to the appropriate regimes: 

x < xw region A , (269) 

x > xw region B . (270) 



The energy distribution of polarized r lepton reads 



dx 



= 12V 



\f(y)±Af(x) 



(271) 
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The function /(x) represents the unpolarizcd case, 

/(*) = ^V* 2 - 4?7 {C [x 2 - 3y(l + 8 V ] + (3x - 6tj)(2 - ar)} , (272) 
with 

C = l-T— ■ (273) 

1 + ?7 — x 



The function A/(x) reads, 



1 



with 



A/(x) = -y/r,\x - 1\ fa* + -r, fa , (274) 



-5A 3 /(x - l) 4 + 3A(4t ? — A - A 2 )/(a; - l) 3 + (4t?A - 4t? + A 
+7A 2 + A 3 )/(x - l) 2 + (-1 + 4r/A - 28// + 15A - A 2 - A 3 )/(x - 1) 
-1 + I2x<q - llxX + 7x - x 2 + Ur/X - 24?y + 14A - 11A 2 , (275) 



j A 



y/x 2 - 4?7 [ 15A 2 £/(x - l) 3 + (IOt^ 2 - 16r?£ + 24A£ - 10A£ 2 
-20A - 6\ 2 £)/(x - l) 2 + (-4 - 14r?A£ 2 - 48t/£ + 66r^ 2 - 24r?£ 3 
+8?y 2 £ 3 - 76A£ + 14A£ 2 + 48A + 3A 2 £ + 25^ - 26£ 2 + 8£ 3 )/(x - 1) 
+3 - 3a; + 57r?£ - 22?y£ 2 - 12A£ - 21A + 34£ - 18£ 2 + 8£ 3 ] , (276) 

cj)B = 15CA 2 /(a;-l) 3 + (60?7CA-1677C-3077C 2 A-16CA-21CA 2 

+ 10C 2 A)/(.t - l) 2 + (-104?7CA - 8477C + 52ryC 2 A + 122r;C 2 - 40?7C 3 
+ 160?7 2 C - 160r; 2 C 2 + 40?7 2 C 3 - 24(A + 9CA 2 + 17( - 4( 2 A - 22( 2 
+8C 3 )/(x - 1) + 18 - 29xr/ + 27xA - 21x + 3x 2 + 46?]CA - 59r/C 
-14?7C 2 A + 78ryC 2 - I677C 3 - 71ryA + 59?? - 43?7 2 C + 26?? 2 C 2 - 8r/ 2 C 3 
+46?7 2 - 6CA - 3CA 2 + 69( - 6C 2 A - 52C 2 + 16C 3 - 42A + 24A 2 (277) 

where 

£ = 2-C, X = P + V ■ (278) 
The function * can be written in the form, 

_ f arccosw mm - arccosw mQ:!; , x < 1 , , 

\ arcoshw maa; - arcosha; mm , x > 1 

with 

_ 2(x l)^mm,maa: ~t~ x(x m x) 2fj /OQfA 

x\j(x m - x) 1 + \r\p 



Due to terms containing inverse powers of [x — 1) the expression fl274| ) for A/(x) 
is apparently divergent at x = 1. However, expanding A/(x) in powers of (x— 1) 
for x < 1 and x > 1 one can check that this function is regular at x = 1. 
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Figure 21: Integrated polarisation of r lepton along the direction of virtual W 
in the Born approximation (dot-dashed), including HQET corrections (long- 
dashed), perturbative corrections (short-dashed) and with both corrections in- 
cluded (solid) as functions of the scaled r energy x. The mass of the b quark 
taken to be 4.75 GeV, c quark 1.35 GeV. The strong coupling constant a s — 0.2. 
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Figure 22: Polarisation of r lepton along the direction of virtual W in the Born 
approximation (dashed) and including one loop corrections (solid) as functions 
of the scaled r energy x. The mass of the b quark taken to be 4.75 GeV, c quark 
1.35 GeV. a s = 0.2. 



The HQET contribution to the decay distributions is known to render them 
unreliable near the endpoint values of the tauon energy. This ambiguity reveals 
itself in the polarisation as well. Similar problems appear also in calculations of 



perturbative corrections! 32 5q] . All these problems are cured if instead of dis- 
tributions their moments are considered |(3(], |9[ [35j. In the case of r polarisation 
a better defined quantity is the integrated polarisation 

p ^w= dx \^-^) dx \^ x - + ^) (281) 

where both the lowest-order perturbative and the HQET terms are included. In 
Fig.^l] the integrated polarisation is shown as a function of the scaled energy x 
of the r lepton. The lowest order prediction corresponds to the dot-dashed line 
and the solid line is obtained including HQET and one-loop perturbative cor- 
rections. The latter have been evaluated numerically. The diagram also shows 
the effect of the perturbative corrections only. The perturbative correction to 
the polarisation is regular so there is no need to integrate over the scaled en- 
ergy x. Fig. [2^ shows the differential polarisation at tree level and including 
the one-loop correction. It is remarkable that this correction, though smaller 
in weight than the correction to the decay rate, is not suppressed to the extent 
that it is for the longitudinal polarisation. On the other hand, we see the two 
polarisations coincide at the endpoint where the energy of the tauon is large so 
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that its direction is almost equal to that of the pair of leptons. This is of course 
expected. 

On integration over the whole range of the charged lepton energy one arrives 
at the total polarisation at the tree level corrected for the 0(1/ mf) effects as 
predicted by HQET. For m b = 4.75 GeV and m c = 1.35 GeV, we obtain 

P = -0.7235 + 4.21— \ + 1.48— % . (282) 

m b m b 

Taking K b = 0.15 GeV 2 , G b = -0.18 GeV 2 and a s = 0.2 we obtain P = -0.709. 

Although we are mostly concerned with the tau lepton polarisation here, 
the formulae derived in the present work may well be used in evaluating the 
polarisation of the light leptons. Interestingly, in the limit of a vanishing mass 
of the cha rged l epton the polarisation falls to zero apart from the endpoints, 



c.f. ([274]) and ( 253 ). It is due to th e ch iral V — A structure of the weak 
charged current that, according to Eq. ( |l45[ ) , the decay widths with a definite 
polarisation differ by a term proportional to m T s tl . The polarisation four- vector 
of the charged lepton can be decomposed as follows: 

E 



= (s°,s) = iVl-(s-L)V-L,-Vl-(s±) 2 , (283) 



777 777 



where s± is understood to mean the part of the three- vector s perpendicular to 
the direction of the charged lepton. The quantities E and p denote, respectively, 
the energy and the three-momentum value of the charged lepton. This form can 
easily be seen to meet the definition of the polarisation four-vector. As the 



lepton mass approaches zero Eq.(283) gives 



ms^ PS y/i - (Sjj 2 + 777 (0, Sj_, 0) . (284) 

However, if we want to keep the angle subtended by the polarisation vector and 
the lepton momentum constant the parallel part of the polarisation should be 
proportional to the perpendicular one, thereby forcing the factor of yT- (s±) 2 



to be of order of m/E. Then the r.h.s of Eq.(284) tends to zero for m — > 0. 
For the vanishing charged lepton mass the polarisation can be non-zero only 
where the virtual W boson is collinear with the charged lepton. In general the 
contribution to polarisation is appreciable only for W direction within the cone 
defined by the condition 

\s x \/\s\=0(m/E) . (285) 
In particular this happens if p is much larger than the energy of the neutrino. 



For semitauonic B decays the condition (285) is satisfied in the whole phase 



space and the resulting polarisation is fairly large. 
7.6 Summary 

We have calculated the Born term and HQET first power correction to the 
polarisation of the r meson with respect to the momentum of the intermediate 
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W boson. They are presented in an analytical form. Apart from that, the 
one-loop perturbative corrections have been evaluated numerically. In contrast 
to the longitudinal one, this polarisation receives noticeable corrections from 
perturbation expansion. All of the considered contributions have been presented 
in a diagram. 
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8 Summary 



We have calculated a few quantities characterising the semileptonic decay of 
the B mesons. Perturbative and nonperturbative corrections have been found 
to the double distribution in terms of hadronic mass and electron energy for the 
decay B — > X u and to the polarisation of the r lepton in the decay B — > X c . 

The longitudinal polarisation of the r lepton receives only a tiny correction 
at one-loop level and it is nearly equal to the Born result to an excellent ap- 
proximation. As such, this quantity may be of use in determining the quark 
mass differences. 

The polarisation of r is best measured with respect to the momentum of the 
intermediate W boson, so we have studied the appropriate distributions. We 
have found the leading HQET corrections of order 1 / m| analytically as well as 
computed numerically the one- loop perturbative corrections. The latter have 
turned out not to be negligible, contrary to the situation with the longitudinal 
polarisation. 

The problem of measuring the V u b matrix element has been addressed. We 
have proposed a quantity to be compared with experimental data and found the 
theoretical prediction for it including one-loop perturbative and leading twist 
nonperturbative corrections. The estimated error on the extracted value is 10%. 

The results presented in this paper, save for the HQET corrections to the 
longitudinal polarisation, have been found for the first time in the papers 149]- 
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